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Abstract 

Koide's mass formula is an empirical relation among the charged lepton masses which holds 
with a striking precision. We present a model of charged lepton sector within an effective field 
theory with U(3) x SU(2) family gauge symmetry, which predicts Koide's formula within 
the present experimental accuracy. Radiative corrections as well as other corrections to 
Koide's mass formula have been taken into account. We adopt a known mechanism, through 
which the charged lepton spectrum is determined by the vacuum expectation value of a 
9-component scalar field On the basis of this mechanism, we implement the following 
mechanisms into our model: (1) The radiative correction induced by family gauge interaction 
cancels the QED radiative correction to Koide's mass formula, assuming a scenario in which 
the U(3) family gauge symmetry and SU(2)i weak gauge symmetry are unified at 10 2 - 
10 3 TeV scale; (2) A simple potential of $ invariant under U (3) x SU (2) leads to a realistic 
charged lepton spectrum, consistent with the experimental values, assuming that Koide's 
formula is protected; (3) Koide's formula is stabilized by embedding U(3) x SU(2) symmetry 
in a larger symmetry group. Formally fine tuning of parameters in the model is circumvented 
(apart from two exceptions) by appropriately connecting the charged lepton spectrum to 
the boundary (initial) conditions of the model at the cut-off scale. We also disucss some 
phenomenological implications. 



1 Introduction 



Among various properties of elementary particles, the spectra of the quarks and leptons exhibit 
unique patterns, and their origin still remains as a profound mystery. Within the Standard Model 
(SM) of elementary particles, the origin of the masses and mixings of the quarks and leptons is 
attributed to their interactions with the (as yet hypothetical) Higgs boson. Namely, these are the 
Yukawa interaction in the case of the charged leptons and quarks, and possibly the interaction 
represented by dimension-5 operators in the case of the left-handed neutrinos. Even if these 
interactions will be confirmed experimentally in the future, since the couping constants of these 
interactions are free parameters of the theory, the underlying mechanism how the texture of these 
couplings is determined would remain unrevealed. 

There have been many attempts to approach the mystery of the fermion masses by identifying 
empirical relations among the observed fermion masses and exploring underlying physics that 
would lead to such relations. In particular, Koide's mass formula is an empirical relation among 
the charged lepton masses given by [1] 



which holds with a striking precision. In fact, substituting the present experimental values 
of the charged lepton masses [2], the formula is valid within the present experimental accura- 
cies. The relative experimental error of the left-hand side (LHS) of eq. (1) is dominated by 
(m ji /m T ) 1 '' 2 (Am T /m T ) (Am T is the experimental error of m T ) and is of order 1CT 5 . A simple 
mnemonic of the relation (1) is that the angle between the two vectors {^frn~ e , ^ m ni a/"v) an d 




Given the remarkable accuracy with which Koide's mass formula holds, many speculations have 
been raised as to existence of some physical origin behind this mass formula [4, 3, 5, 6, 7, 8, 9]. 
Despite the attempts to find its origin, so far no realistic model or mechanism has been found 
which predicts Koide's mass formula within the required accuracy. The most serious problem one 
faces in finding a realistic model or mechanism is caused by the QED radiative correction [7]. 
Even if one postulates some mechanism at a high energy scale that leads to this mass relation, 
the charged lepton masses receive the 1-loop QED radiative corrections given by 



fh(fi) and m pole denote the running mass defined in the modified-minimal-subtraction scheme (MS 
scheme) and the pole mass, respectively; /i represents the renormalization scale. It is the pole mass 
that is measured in experiments. Suppose fhi(fi) (or the corresponding Yukawa couplings ^(/x) ) 
satisfy the relation (1) at a high energy scale fi 3> M-^. Then mf olc do not satisfy the same relation 
[6, 7]: Eq. (1) is corrected by approximately 0.1%, which is 120 times larger than the present 
experimental error. Note that this correction originates only from the term — 3a /(4n) x fhi log(mf ) 
of eq. (2), since the other terms, which are of the form const, x fhi, do not affect the relation (1). 
This is because, the latter corrections only change the length of the vector (^/m^, ^/m^, v /m T ) 
but not the direction. We also note that log(m^) results from the fact that fhi plays a role of an 
infrared (IR) cut-off in the loop integral. 

The 1-loop weak correction is of the form const, x fhi in the leading order of fhj/M^ expan- 
sion; the leading non-trivial correction is (9(GVmf/7r) whose effect is smaller than the current 




(1) 




(1,1,1) equals 45° [3]. 
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experimental accuracy. Other radiative corrections within the SM (due to Higgs and would-be 
Nambu-Goldstone bosons) are also negligible. 

Thus, if there is indeed a physical origin to Koide's mass formula at a high energy scale, we 
need to account for a correction to the relation (1) that cancels the QED correction. Since such 
a correction is absent up to the scale of O(Mw) to our present knowledge, it must originate from 
a higher scale. Then, there is a difficulty in explaining why the size of such a correction should 
coincide accurately with the size of the QED correction which arises from much lower scales. There 
are also other less serious, but important questions that are often asked: (1) Why do not quark 
masses satisfy the same or a similar relation? (2) In Koide's formula the three lepton masses 
appear symmetrically. Then why is there a hierarchy among these masses, m e <C m M <C m r ? If 
there is indeed a physical origin to Koide's mass formula, there must be reasonable answers to all 
of these questions. 

Among various existing models which attempt to explain origins of Koide's mass formula, we 
find a class of models particularly attractive [10, 13]. These are the models which predict the mass 
matrix of the charged leptons to be proportional to the square of the vacuum expectation value 
(VEV) of a 9-component scalar field (we denote it as $) written in a 3-by-3 matrix form: 

M e cx ($>($>. (3) 

Thus, (y/m e , y/rn^, ^m T ) is proportional to the diagonal elements of (<£>) in the basis where it is 
diagonal. The VEV ($) is determined by minimizing the potential of scalar fields in each model. 
Hence, the origin of Koide's formula is attributed to the specific form of the potential which 
realizes this relation in the vacuum configuration. Up to now, no model is complete with respect 
to symmetry: Every model requires either absence or strong suppression of some of the terms in 
the potential (which are allowed by the symmetry of that model), without justification. 

In this paper, we study possible connections between family (horizontal) gauge symmetries 
and Koide's formula and the charged lepton spectrum. These will be discussed within the context 
of an effective field theory (EFT) which is valid below some cut-off scale. In particular we address 
the following points: 

(i) We propose a possible mechanism for cancellation of the QED radiative correction to Koide's 
mass formula. 

(ii) We propose a mechanism that produces the charged lepton spectrum, which is hierarchical 
and approximates the experimental values, under the assumption that Koide's formula is 
protected by some other mechanism. 

(iii) We present a model of charged lepton sector based on £7(3) x SU(2) family gauge symmetry, 
incorporating the mechanisms (i)(ii). A new mechanism that stabilizes Koide's formula is 
incorporated in this model. 

(Among these, we have reported the main point of (i) separately in [11].) 

In our study we adopt the mechanism eq. (3) for generating the charged lepton masses at tree 
level of EFT, for the following reasons. First, the mechanism allows for transparent and concise 
perturbative analyses of models, which is crucial in keeping radiative corrections under control. 
This may be contrasted with models with other mass generation mechanisms, such as dynami- 
cal symmetry breaking or composite lepton models, which typically involve strong interactions. 
Secondly, since $ is renormalized multiplicatively, the structure of radiative corrections becomes 
simple, as opposed to cases in which VEVs of more than one scalar fields contribute to the charged 
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lepton spectrum. In short, this type of mass generation mechanism is pertinent to serious analyses 
of radiative corrections to Koide's formula, which is a distinguished aspect of this study. 

We alert in advance that we do not solve the hierarchy problem or fine tuning problem of the 
electroweak scale. We cannot explain how to stabilize the electroweak symmetry-breaking scale 
against other higher scales included in our model. Solution to this problem is beyond the scope 
of this paper. 

The paper is organized as follows. In Sec. 2, we explain philosophy of our analysis using EFT 
and argue for its validity and usefulness. We also give a brief overview of the ideas presented in 
this paper. In Sec. 3, we explain the mechanism for cancelling the QED corrections to Koide's 
formula. In Sec. 4, we present a potential for generating a realistic charged lepton spectrum, 
assuming that Koide's formula is protected. In Sec. 5, we analyze a minimal potential whose 
vacuum corresponds to a desired lepton spectrum. In Sec. 6, we extend the potential by including 
another field. In Sec. 7, we introduce a higher-dimensional operator which generates the lepton 
masses and compute corrections to Koide's formula. In Sec. 8, we discuss the energy scales and 
unsolved questions in our model. In Sec. 9, we discuss phenomenological implications of our model. 
In Sec. 10 summary and discussion are given. Technical details are collected in Appendices. 



2 EFT Approach and Brief Overview of the Model 

Throughout this paper, we consider an EFT which is valid up to some cut-off scale denoted by 
A Mw)- In this EFT, we assume that the charged lepton masses are induced by a higher- 
dimensional operator 

O = ipu ®ik ®kj (p e Rj (4) 

(or by other similar operators, as will be described later). Here, ipLi = {yhii^Li) T denotes the 
left-handed lepton SU (2) l doublet of the i-th generation; e^j denotes the right-handed charged 
lepton of the j-th generation; ip denotes the Higgs doublet field. They are respectively assigned 
to the standard representations of the SM gauge group. By contrast, a 9-component scalar field 
$ is absent in the SM and a singlet under the SM gauge group. We suppressed all the indices 
except for the generation (family) indices i,j,k = 1,2,3. (Summation over repeated indices is 
understood throughout the paper unless otherwise stated.) The dimensionless Wilson coefficient 
of this operator is denoted as «(//). Once $ acquires a VEV, the operator O will effectively 
be rendered to the Yukawa interactions of the SM; after the Higgs field also acquires a VEV, 
((f) = (0, v ew / V2) T with v ew 250 GeV, the operator will induce the charged-lepton mass matrix 
of the form eq. (3) at tree level: 

^ ree ~(*>($>. (5) 

For a moment, let us assume that the dimension-4 Yukawa interactions ^uWRj are prohibited 
by some mechanism. This will be imposed explicitly by a symmetry in our model to be discussed 
through Sees. 3-9. 

We now explain philosophy of our analysis using EFT. Conventionally a more standard ap- 
proach for explaining Koide's mass formula has been to construct models within renormalizable 
theories. Nevertheless, the long history since the discovery of Koide's formula shows that it is 
quite difficult to construct a viable renormalizable model for explaining Koide's relation. It is 
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Figure 1: Diagram which induces the higher-dimensional operator O = ^pLi &ik &kj peRj- 

likely that we are missing some essential hints to achieve this goal, if the relation is not a sheer co- 
incidence. In this paper we will show that, within EFT, explanation of Koide's formula is possible 
by largely avoiding fine tuning of parameters. Consistency conditions (with respect to symmetries 
of the theory) can be satisfied relatively easily in EFT, or in other words, they can be replaced 
by reasonable boundary conditions of EFT at the cut-off scale A without conflicting symmetry 
requirements of the theory. (See Sec. 5.) Even under this less restrictive theoretical constraints, 
we may learn some important hints concerning the relation between the lepton spectrum and 
family symmetries. These are the role of specific family gauge symmetry in canceling the QED 
correction, the role of family symmetry in stabilizing Koide's mass relation, or the role of family 
symmetry in realizing a realistic charged lepton spectrum consistently with experimental values. 
These properties do not come about separately but are closely tied with each other. These features 
do not seem to depend on details of more fundamental theory above the cut-off scale A but rather 
on some general aspects of family symmetries and their breaking patterns. Thus, we consider that 
our approach based on EFT would be useful even in the case in which physics above the scale A 
is obscure and may involve some totally unexpected ingredients - - as it was the case with chiral 
perturbation theory before the discovery of QCD. 

Before discussing radiative corrections within EFT, one would be worried about effects of 
higher-dimensional operators suppressed in higher powers of 1/A. Indeed, using the values of tau 
mass and the electroweak symmetry breaking scale t> ew , one readily finds that v^/A ^ 0.1 (vi are 
the diagonal elements of ($) in the basis where it is diagonal). Hence, naive dimensional analysis 
indicates that there would be corrections to Koide's formula of order 10% even at tree level. We 
now argue that this is not necessarily the case within the scenario under consideration. We may 
divide the corrections into two parts. These are (i) 1/A n corrections to the operator O of eq. (4) 
(the operator which reduces to the SM Yukawa interactions after $ is replaced by its VEV), and 
(ii) 1/A" corrections to the VEV of $. 

Concerning the corrections (i), we may consider the following example.* Suppose that the 
operator O is induced from the interactions 

£ = y x r, ( «I», ; // /0 + M H Ri H Li + y 2 //,,<I\, //;, ; + M'H' Rl H' Li + y 3 H' LiV e m + (h.c.) (6) 

through the diagram shown in Fig. 1, after fermions Hl,r and H' L R have been integrated out. 
Fermions H^ ji and H' L R are assigned to appropriate representations of the SM gauge group such 
that the above interactions become gauge singlet. For instance, in the case that v 3 /M' > 3, 
2/1,2/2, 2/3 ~ 1 and v cw /M' < 3 x 10~ 3 , one finds, by computing the mass eigenvalues, * that the 

*A similar mechanism is used in [10, 12] to induce O; corrections by higher-order terms in 1/A™ have not been 
discussed, however. 

' Since the values of m T and v cw are known, once we choose the values of vz/M'( > 3) and yi,2/2,2/3(~ 1), the 
value of vs/M( < 0.03) will be fixed. Then the mass eigenvalues corresponding to the SM charged leptons can be 
computed in series expansion in the small parameters v cw /M', Vi/M and vf/(MM') = \/2mi/v cw . 
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Figure 2: EFT 1-loop diagrams which generate higher-dimensional operators contributing to the charged 
lepton spectrum. Dashed line represents <J?; (8> represents the higher-dimensional operator which generate 
the charged lepton masses at tree-level [corresponding to O of eq. (4)]. 

3 3 

largest correction to the lepton spectrum eq. (5) arises from the operator ^l^ 6( P^r'j its 

_ 1 VJ. 1 VJ. 

contribution to the tau mass is 5m T /m T = (m T /v cw ) 2 ~ 5 x 1CT 5 . This translates to a correction 
to Koide's relation of 3 x 1CT 6 , due to the suppression factor 4g (m At /m r ) 1 / 2 (<5m T /m T ). Thus, 
this is an example of underlying mechanism that generates the operator O without generating 
higher-dimensional operators conflicting the current experimental bound. If we introduce even 
more (non-SM) fermions to generate the leading-order operator O, one can always find a pattern 
of spectrum of these fermions, for which higher-dimensional operators are sufficiently suppressed, 
since the number of adjustable parameters increases. (Another example of underlying mechanism 
may be the one proposed in [8], based on the idea of [14].) In general, sizes of higher-dimensional 
operators depend heavily on underlying dynamics above the cut-off scale. 

Let us restrict ourselves within EFT. If we introduce only the operator O, by definition this 
is the only contribution to the charged lepton spectrum at tree level. Whether loop diagrams 
induce higher-dimensional operators which violate Koide's relation is an important question, and 
a detailed analysis is necessary. This is the subject of the present study, where the result depends 
on the mechanisms how Koide's formula is satisfied and how the charged lepton spectrum is 
determined, even within EFT. The conclusion is as follows. Within the model to be discussed 
in Sees. 3-8, the class of 1-loop diagrams shown in Fig. 2 do not generate operators that violate 
Koide's relation sizably; see Sec. 7. (There is another type of 1-loop diagrams that possibly cancels 
the QED correction; see Sec. 3.) In fact, we do not find any loop-induced higher-dimensional 
operators, which violate Koide's relation in conflict with the current experimental bound. 

Concerning the corrections (ii), in our analysis we introduce specific family gauge symmetries 
and their breaking patterns such that the corrections (ii) are suppressed. 

Since the above example of underlying mechanism that suppresses higher-dimensional oper- 
ators is simple, and since suppression of loop-induced 1/A n corrections within EFT provides a 
non-trivial cross check of theoretical consistency, we believe that our approach based on EFT has 
a certain justification and would be useful as a basis for considering more fundamental models. 

In the rest of this section, we present a brief overview of the basic ideas of the analysis to be 
given through Sees. 3-8, in order to facilitate reading. Our analysis starts from investigating a 
possibility that the radiative correction generated by a family gauge symmetry cancels the QED 
correction to Koide's formula (Sec. 3). We find that U(3) ~ SU(3) x U(l) family gauge symmetry 
has a unique property in this regard. In fact, if ipL and e# are assigned to mutually conjugate 
representations of this symmetry group, the U(3) radiative correction has the same form as the 
QED correction with opposite sign. In particular, if the gauge coupling of U(3) family symmetry 
a F = 5 , f/(47t) satisfies the relation a(m T ) ~ \(Xf(9fV3), both corrections cancel. We speculate 
that this relation would be realized within a scenario in which U(3) family gauge symmetry is 
unified with SU(2)l gauge symmetry at 10 2 -10 3 TeV scale, although we need to fine tune the 
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unification scale within an accuracy of factor 3. 

The non-trivial form of the radiative correction by the U(3) gauge interaction is dictated by the 
£7(3) symmetry and its breaking pattern induced by the VEV (<&). In particular, multiplicative 
renormalizability of ($) ensures that the correction to Koide's formula is independent of the 
renormalization scale /i of the effective potential of $. Namely, the charged lepton pole masses 
are determined, up to a common multiplicative constant, directly by the form of the effective 
potential renormalized at an arbitrary high scale \x (< A), and we may ignore the QED and U(3) 
radiative corrections altogether. For our purpose, it is most convenient to take this scale to be 
fi = A. In this part of our analysis, we assume that (<£>) can be brought to a diagonal form by 
symmetry transformation, and also that Koide's relation for the diagonal elements, 



In the second step, we search for an effective potential for which the eigenvalues of (<&) satisfy 
the relation (7) and reproduce the experimental values of the mass ratios v± : V2 V3 — ^frn~ e : 




^/m^ : yJm T (Sees. 4 and 5). If we choose the renormalization scale to be /i — A, radiative 
corrections to the effective potential essentially vanish within EFT, or in other words, the form 
of the effective potential at this scale is determined by physics above the scale A as boundary 
(initial) conditions of EFT. Hence, our goal is to find an effective potential which satisfies the 
boundary conditions without conflicting symmetry requirements of the theory. Although it may 
seem an easy task, it still involves fairly non-trivial analyses. 

We impose U(3) x SU(2) family symmetry as a symmetry of EFT. The motivation of this 
choice is that it is the symmetry possessed by the simplest higher-dimensional operator analyzed 
in the first step. It turns out, however, that this symmetry is not large enough to constrain the 
form of the effective potential sufficiently We therefore further assume a symmetry enhancement. 
Namely, we assume that above the cut-off scale A there is an SU(9) x (7(1) gauge symmetry, and 
this symmetry is spontaneously broken to U (3) x SU(2) below the cut-off scale. The symmetry 
SU(9) x U(l) is motivated by a geometrical interpretation of Koide's relation eq. (7). Within 
this scenario, we still need to introduce an additional scalar field X in order to realize a desirable 
vacuum configuration. Thus, we analyze the vacuum of the general potential of $ and X . (Details 
of the analysis are rather technical.) The conclusion is that in a finite region of the parameter space 
of the potential, Koide's relation is satisfied by the eigenvalues of (<&). Furthermore, the eigenvalues 
can be made consistent with the experimental values of the charged lepton masses without fine 
tuning of parameters. These are realized in the case that certain hierarchical relations among the 
parameters of the potential are satisfied, and these relations do not conflict the requirement of the 
assumed symmetry and symmetry enhancement. We speculate on possible physics scenario above 
the cut-off scale that may lead to (part of) these hierarchical relations. 

So far, these desirable features are satisfied by the eigenvalues of ($). There remains, however, 
a problem that (<&) cannot be brought to a diagonal form by the U (3) x SU{2) symmetry transfor- 
mation, and this contradicts the assumption made in the first step. To remedy this difficulty, we 
introduce yet another field £y such that, with an appropriate potential with $, it can generate an 
appropriate higher-dimensional operator necessary to produce the charged lepton masses. (Sees. 6 
and 7.) Although the potential and the higher-dimensional operator involving £y do not conflict 
the requirement of the assumed symmetry and symmetry enhancement, these would be the most 
unsatisfactory part of our model. This is because it is difficult to speculate any plausible scenario 
above the cut-off scale, which would lead to these potential and higher-dimensional operator. 




(7) 



is satisfied. 
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Figure 3: Diagram for the Hoop correction by the family gauge bosons to the operator O when both 
■ipL and en are in the 3 of 577(3). The diagram on the right-hand-side shows the flow of family charge 
in the leading contribution of the 1/Np expansion (Np = 3); closed loop corresponds to tr(($)($)). 

With all these setups, it is possible to compute the radiative corrections which are induced 
by the diagrams shown in Fig. 2. Due to the specific form of the effective potential of $ and 
X, corrections to Koide's formula turn out to be quite suppressed, as long as the aforementioned 
hierarchical conditions between parameters of the potential are satisfied. As already mentioned, 
this serves as a non-trivial consistency check of the model as an EFT. 

There are a few unsolved questions and incompleteness of the present model and these are 
discussed in Sees. 8 and 10. 

3 Radiative Correction by Family Gauge Interaction 

In this section we introduce family gauge symmetries and consider radiative corrections to the mass 
matrix eq. (5) by the family gauge interaction. First we consider the case, in which the family 
gauge group is 577(3) and both ipp and e# are assigned to the 3 (fundamental representation) of 
this symmetry group. We readily see, however, that with this choice of representation, Koide's 
formula will receive a severe radiative correction unless the family gauge interaction is strongly 
suppressed. In fact, the 1-loop diagram shown in Fig. 3 induces an effective operator 



hence corrections universal to all the charged-lepton masses, (8m e , dm^, 5m T ) oc (1,1,1), are 
induced. This is due to the fact that the dimension-4 operator ipn (p em is not prohibited by sym- 
metry. Here, ap = gp/(4:ir) denotes the gauge coupling constant of the family gauge interaction. 
As noted above, corrections which are proportional to individual masses do not affect Koide's for- 
mula; oppositely, the universal correction violates Koide's formula rather strongly. In order that 
the correction to Koide's formula cancel the QED correction, a naive estimate shows that a^/V 
should be order 10~ 5 , provided that the cut-off A is not too large and that the above operator 
O' is absent at tree level. If O' exists at tree level, there should be a fine tuning between the 
tree-level and 1-loop contributions. The situation is similar if the family symmetry is 0(3) and 
both ip L and e# are in the 3, which is also a typical assignment in existing models. In these cases 
[15] we were unable to find any sensible reasoning for the cancellation between the QED correction 
and the correction induced by family gauge interaction, other than to regard the cancellation as 
just a pure coincidence. Hence, we will not investigate these choices of representation further. 

In the case that tpL is assigned to 3 and e R to 3 (or vice versa) of U(3) family gauge group, 
(i) the dimension-4 operator 4>Li¥ e Ri is prohibited by symmetry, and hence corrections universal 
to all the three masses do not appear; and (ii) marked resemblance of the radiative correction to 
the QED correction follows. We show these points explicitly in a specific setup. 




A 2 



(8) 
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We denote the generators for the fundamental representation of £7(3) by T a (0 < a < 8), 
which satisfy 

tr (T a T p ) = -5 al3 ; T a = T a] . (9) 

T° is the generator of U(l), hence it is proportional to the identity matrix, while T a (1 < a < 8) 
are the generators of SU (3). Here and hereafter, a, f3, 7, . . . represent C/ (3) indices 0, . . . , 8, while 
a,b,c, . . . represent SU(3) indices 1, . . . , 8. The explicit forms of T a are given in Appendix A. 

We assign ipi to the representation (3,1), where 3 stands for the SU(3) representation and 
1 for the U(l) charge, while e R is assigned to (3,-1). Under U(3), the 9-component field $ 
transforms as three (3, l)'s. Explicitly the transformations of these fields are given by 

^l^U^l, e R ^U*e R , ; U = exp (i9 a T a ) , UU* = 1. (10) 

We assume that the charged-lepton mass matrix is induced by a higher-dimensional operator 
similar to O in eq. (4). We further assume that (<£>) can be brought to a diagonal form in an 
appropriate basis. Thus, in this basis after $ and tp acquire VEVs, turns to the lepton mass 
terms as 



m 



tree 

'e 











m tree 











m tree 



OW^i> L MT e e R , Mf ee = I I =^^$ d (/i) 2 , (11) 

where 

$ d (/i)= I v 2 (n) I , ^(/i)>0. (12) 




When all Vi are different, U(3) symmetry is completely broken by (<&) = $d, and the spectrum of 
the U(3) gauge bosons is determined by $ d . 

Note that the operator O in eq. (4) is not invariant under the Z7" (3) transformations eq. (10). 
As an example of one may consider 

oFU^fcMTye,. (13) 

It is invariant under a larger symmetry U(3) x SU(2), under which $ transforms as $ — > U&0 T 
(00 T = 1). In this case, we need to assume, for instance, that the SU(2) symmetry is gauged 
and spontaneously broken at a high energy scale before the breakdown of the U(3) symmetry, 
in order to eliminate massless Nambu-Goldstone bosons and to suppress mixing of the U (3) and 
SU(2) gauge bosons. A more elaborate example of the higher-dimensional operator, which is 
consistent with the symmetry and satisfies eqs. (11) and (12), will be given in Sec. 7. In any 
case, the properties of given by eqs. (11) and (12) are sufficient for computing the radiative 
correction by the U(3) gauge bosons to the mass matrix, without an explicit form of O^. 
We take the U(l) and SU(3) gauge coupling constants to be the same: 

Oiu(i) = a>su(3) = oc F . (14) 
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eR 

(a) 






eR 



eR 



(c) 



Figure 4: 1-loop diagrams contributing to 5m\ oe when ipi and are in the (3, 1) and (3, —1), respec- 
tively, of SU(3) x U(l). The right-hand-sides show flows of family charge, (a) Correction of the form 
ipL eR'- charge flow is connected in one line, showing multiplicative renormalization, (b) correction of 
the form -tpL ^ Z^ipL, and (c) correction of the form eR$Z e eR. 



We compute the radiative correction in Landau gauge, which is known to be convenient for 
computations in theories with spontaneous symmetry breaking. From the diagrams shown in 
Figs. 4(a) (b)(c), we find 

mi(fi), (15) 

(16) 

Here, c is a constant independent of i. The Wilson coefficient k^\^) is defined in MS scheme. 
fi(/i)'s are defined as follows: The VEV of $ at renormalization scale //, = ($(//)) given by 

eq. (12), is determined by minimizing the 1-loop effective potential in Landau gauge (although we 
do not discuss the explicit form of the effective potential in this section); $ is renormalized in MS 
scheme. We ignored terms suppressed by m 2 /v 2 (<^i 1) in the above expression. Note that the pole 
mass is renormalizat ion-group invariant and gauge independent. Therefore, the above expression 
is rendered gauge-independent if we express Vi(fj) in terms of gauge-independent parameters, such 
as coupling constants defined in on-shell scheme. 

The coefficient of log fi 2 is determined by the sum of the anomalous dimension of the Wilson 
coefficient and twice of the wave-function renormalization of <£>. (The former is gauge 

independent, while the latter is not.) The term \ogvf originates from the role of the gauge boson 



r pole 

dm- = 



3 ap 
'~8tT 



log 



rriiilj) = 



K 



(n)v e 



V2A 2 



Vi{fl) 2 . 



+ c 
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masses as an IR cut-off of the loop integral, hence it reflects the spectrum of the gauge bosons. 
The sign in front of log fi 2 is opposite to that of the QED correction eq. (2), which results from the 
fact that ipL and have the same QED charges but mutually conjugate (opposite) Z7" (3) charges. 

In Landau gauge, the diagrams in Figs. 4(b) (c) are finite and flavor independent, i.e. propor- 
tional to 5ij in terms of the family indices; hence they contribute only to the constant c. Apart 
from this constant, the difference between the QED correction and the correction (15) resides in 
the factors 



a rrii log ( ) $ij 



vs. 



(17) 



in the QED self-energy diagram and the diagram in Fig. 4(a), respectively. (No sum over % is 
taken in the former factor.) One may easily identify the factor 2 difference in the coefficients of 
log ji 2 using the Fierz identity 



Ik 



(18) 



From this identity, it follows that the operator is multiplicatively renormalized; see family 
charge flow in Fig. 4(a). (Mp) a p denotes the mass matrix of the family gauge bosons. After 
diagonalization, one obtains the spectrum of family gauge bosons as 



gl 



i (^) 2 + vi + - 2 {v\ + vim^r + (^ 4 ) 2 } + vi (^) 



7 b\2 



+ kvl + Vl){(^) 2 + (^ ? ) 2 } + + -3){(^ 8 ) 2 + (^) 2 ) 



• (19) 



The mass eigenstates JF* are labelled in the order of their masses, which are given by 

f _ ./O f 8 

J fJ, V J fl 



fO 



fi fS fO fi fS 



<p3,4 



f, 



1,2 



-r-6,7 



4,5 



778,9 



f6,7 
J u 



^3 



(20) 



Hence, 



f a T a = 



1 -r-l 



1 -172 _ i I f, 

y/2 V 2 y^t 



78 



1 -r-5 



(21) 



The form of the radiative correction given by eqs. (15) and (16) is constrained by symmetries 
and their breaking patterns. As the diagonal elements of the VEV, 1*3 > v-i > v\ > 0, are 
successively turned on, gauge symmetry is broken according to the pattern: 



U(3) -> U(2) -> U{1) -> nothing. 



(22) 
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At each stage, the gauge bosons corresponding to the broken generators acquire masses and 
decouple. Furthermore, the vacuum Q d and the family gauge interaction respect a global U(l)vi x 
U(l)v2 x U(l)v3 symmetry generated by 

ipL * UdipLi &r > U* d e R , $ d - U d U* d (23) 

with 



1<P\ 








e i4>2 o 





e^ 3 



U d =\ ; 0,GR. (24) 

\ e i<fe / 

The operator after symmetry breakdown, eq. (11), is not invariant under this transformation 
but the variation can be absorbed into a redefinition of v^s. As a result, the lepton mass matrix 
has a following transformation property: 



Me 



= U d M e U* d . (25) 

*Vi exp(i<f>i) 



This is satisfied including the 1-loop radiative correction. The symmetry breaking pattern eq. (22) 
and the above transformation property constrain the form of the radiative correction to Sm^ ole oc 
fj 2 [log(|f j| 2 ) + const.], where the constant is independent of %. Note that |t>j| 2 in the argument of 
logarithm originate from the gauge boson masses, which are invariant under Vi — > Vi exp(i</>j) . 

The universality of the SU(3) and U(l) gauge couplings eq. (14) is necessary to guarantee the 
above symmetry breaking pattern eq. (22). One may worry about validity of the assumption for 
the universality, since the two couplings are renormalized differently in general. The universality 
can be ensured approximately if these two symmetry groups are embedded into a simple group 
down to a scale close to the relevant scale. There are more than one ways to achieve this. A 
simplest way would be to embed SU (3) x U(l) into SU (4). It is easy to verify that the 4 of SU (4) 
decomposes into (3, — \) © (1, |) under SU(3) x U(l). Hence, the 6 (second-rank antisymmetric 
representation) and 6 of SU(A), respectively, include (3, —1) and (3, 1). 

Within the effective theory under consideration, the QED correction to the pole mass is given 
just as in eq. (2) with rhi(n) replaced by rrii(n). Recall that corrections of the form const, x do 
not affect Koide's formula. Then, noting logvf = |logmf + const., one observes that if a relation 
between the QED and family gauge coupling constants 

a = -aF (26) 

is satisfied, the 1-loop radiative correction induced by family gauge interaction cancels the 1-loop 
QED correction to Koide's mass formula. 

In fact, with the relation (26), cancellation holds for all the leading logarithms generated 
by renormalization group: the coefficient of log/z 2 of the QED correction is determined by the 
1-loop anomalous dimension of the running mass, while the coefficient of logyU 2 of eq. (15) is 
determined by the anomalous dimension of the Wilson coefficient and twice of the wave- 
function renormalization of $; they are resummed in the same way by 1-loop renormalization 
group equations. The renormalization group evolution and the symmetry breaking pattern eq. (22) 
in the scale range across the family gauge boson masses dictate how logmf's induced by family 
gauge interaction are resummed. The renormalization group evolution and the same symmetry 
breaking pattern in the QED sector dictate the log mf resummation of the QED correction, in the 
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scale range across the lepton masses. If TOjlogm 2 cancel at 1-loop, log to? dependences in all the 
leading logarithms mi[alog(/j, 2 /m 2 )] n also cancel. On the other hand, effects of the running of a 
and ap do not cancel. It is related to the question which we stated in the Introduction: What are 
the relevant scales for the coupling constants in the relation (26)? The scale of a is determined by 
the lepton masses, while the scale of «f is determined by the family gauge boson masses, which 
should be much higher than the electroweak scale. 
Suppose the relation (26) is satisfied. Then 



holds including the leading logarithms generated by the running of and v^s. This is valid for 
any value of /i. This means, if i>j(//)'s satisfy 



at some scale //, Koide's formula is satisfied at any scale [i. This is a consequence of the fact that 
$ is multiplicatively renormalized. Generally, the form of the effective potential varies with scale 
/i. If the relation (28) is realized at some scale as a consequence of a specific nature of the effective 
potential (in Landau gauge), the same relation holds automatically at any scale. Although these 
statements are formally true, physically one should consider scales only above the family gauge 
boson masses, since decoupling of the gauge bosons is not encoded in MS scheme. For our purpose, 
it is most appropriate to use eq. (27) to relate the charged lepton pole masses with the VEV at 
the cut-off scale, i.e. /i — A, which sets a boundary (initial) condition of the effective theory. 

The advantages of choosing Landau gauge in our computation are two folds: (1) The compu- 
tation of the 1-loop effective potential for the determination of ($) becomes particularly simple 
(as well known in computations of the effective potential in various models); in particular there is 
no 0(cxf) correction to the effective potential. (2) The lepton wave-function renormalization is 
finite; as a consequence, the diagrams in Figs. 4(b) (c) are independent of ($(//)) and independent 
of flavor. Due to the former property, there is no 0(ap) correction to the relation eq. (28) if it is 
satisfied at tree level. Due to the latter property, £mf ole is determined essentially by the diagram 
in Fig. 4(a) and a simple relation to (<&(//)) follows. 

Let us comment on gauge dependence of our prediction. If we take another gauge and express 
the radiative correction 5mf ole in terms of (<&(//)), the coefficient of log(/x 2 / ($) 2 ) changes, and other 
non-trivial flavor dependent corrections are induced. Suppose the relation eq. (28) is satisfied at 
tree level.* The VEV (<&) in another gauge receives an Oiptp) correction, which induces a correction 
to eq. (28) at 0(cxf)- These additional corrections to 5m^ olc at 0(cxf) should cancel altogether 
if they are reexpressed in terms of the tree-level v^s which satisfy eq. (28), since the O(oif) 
correction to the relation (28) vanishes in Landau gauge. General analyses of gauge dependence 
of the effective potential may be found in [17]. 

Now we speculate on a possible scenario how the relation (26) may be satisfied. Since the 
relevant scales involved in a and a? are very different, we are unable to avoid assuming some 
accidental factor (or parameter tuning) to achieve this condition. Instead we seek for an indirect 
evidence which indicates such an accident has occurred in Nature. The relation (26) shows that 
the value of ap is close to that of the weak gauge coupling constant aw, since sin 2 9w(M\y) is close 
to 1/4. In fact, within the SM, ^aw(fj) approximates a(m T ) at scale /i ~ 10 2 -10 3 TeV. Hence, 

*To simplify the argument we consider only those gauges in which tree-level vacuum configuration is gauge 
independent, such as the class of gauges considered in [16]. 



(27) 





(28) 
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if the electroweak £77(2) L gauge group and the Z7(3) family gauge group are unified around this 
scale, naively we expect that 



a « - a F (29) 

is satisfied. Since aw runs relatively slowly in the SM, even if the unification scale is varied within 
a factor of 3, Koide's mass formula is satisfied within the present experimental accuracy. This 
shows the level of parameter tuning required in this scenario. 

We may generalize our setup and see how the radiative correction alters. If ipL and cr are 
assigned to (3,Q^ L ) and (3,Q eR ), respectively, the correction eq. (15) generalizes to 



r pole ®-F 



(Qf L -Qe R + l) log {j^J^j + C ' 



mi(fjt) , (30) 



where d is a flavor-independent constant. Thus, the form rrii logm 2 is maintained. This is not the 
case if we vary the U (1) charge of $, which violates the breaking pattern of gauge symmetry eq. (22) 
strongly.''' The form rriilogm 2 is maintained in yet another generalization, in which U(3) x £7(3) 
symmetry is gauged. We introduce another field X : (1, 0, 6, 2) under SU (3) x U(l) x SU (3) x £7(1). 
The symmetry transformations are given by ipL — * UipL, &r — > U*en, $ — > UQV\ X — > VY<V T 
with U = exp(i9 a T a ), V = exp(i6 a T a ). We assume that (£) = u E l with w s < v 1 ,v 2 ,v 3 , and that 
the lepton masses are generated by a higher-dimensional operator 

0? = ^j L *V* T <pe R . (31) 
For the assignment ipL '■ (3, 1, 1, 0) and e# : (3, — 1, 1, 0), the radiative correction reads 



r pole 

dm- = 



3 at r ' - 2 



8vr a 2 F + a' F 



rriiifi) , (32) 



where ap and a' F denote, respectively, the gauge couplings of the first U(3) and second U(3) 
symmetries. * Thus, the coefficient of \ogm 2 varies in different setups. Accordingly the condition 
for the cancellation of the QED correction changes from eq. (26). We need to seek for other possible 
scenarios which lead to such conditions, or maybe to let the cancellation be a sheer coincidence. 
The level of fine tuning required for the coupling(s) is about 1% to meet the present expermental 
accuracy of Koide's formula. 

In the rest of this paper, we do not consider these generalizations. We adhere to eq. (26), 
assuming the scenario in which SU(2)l and U(3) gauge symmetries are unified at around 10 2 - 
10 3 TeV. In this paper we do not construct a model which incorporates this unification scenario. 
We simply assume that this unification scenario is realized in the underlying full theory, in which 
the unification scale is at or around the cut-off scale A of our effective theory; we further assume 
that the hierarchy between v 3 and A(> v 3 ) is mild; see discussions in Sees. 2 and 8. 



^Even in the case in which only the U(l) charges of ipL and en are varied, this symmetry breaking pattern is 
violated but only softly through the gauge interaction of ipL an d er- By contrast, varying the U(l) charge of $ 
affects the spectrum of the gauge bosons. 

*If ap = a' F , SU(3) x SU(3) x U(1)a can be embedded into SU(6). In this case, tpL and can be assigned 
to the 6 and 6 of SU(6), respectively. The remaining U(l)v, corresponding to the lepton number, is unbroken, so 
it may be taken as a global symmetry. 
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Figure 5: Geometrical interpretation of eq. (34). Eq. (9) defines the inner product in a 9-dimensional 
real vector space spanned by the basis {T a }. Since $°T°, $ a T a and $ = $ Q T a form an isosceles right 
triangle, the angle between T° and $ is 45°. This is Koide's formula in the basis where <J? is diagonal [3]. 



4 Potential Minimum and Charged Lepton Spectrum 

The analysis in the previous section indicates relevance of the U(3) family gauge symmetry in 
relation to the charged lepton spectrum and Koide's mass formula. In this section we study the 
potential of <£> invariant under this family symmetry and its classical vacuum. In particular, we 
propose a mechanism for generating a realistic charged lepton spectrum, assuming that Koide's 
mass relation is protected. For later convenience, we express components of $ using T a , denned 
in eqs. (9) and (98), as the basis: 

$ = $ a T a . (33) 

In general $ a takes a complex value. 

The largest symmetry that can be imposed on the higher-dimensional operator is U (3) x 
£7(3). [An example is given in eq. (31).] We may consider the potential of $ consistent with 
this symmetry, allowing only operators with dimension 4 or less. A general analysis shows that, 
for any choice of the parameters (couplings) of this potential, the classical vacuum ($), after its 
diagonalization, does not satisfy the relation (28) [18]. Namely, there is no classical vacuum that 
leads to Koide's mass formula. If we impose a smaller symmetry on the potential of $, it is 
possible to tune the parameters in the potential and realize the relation (28) as well as a realistic 
charged lepton spectrum. We were, however, unable to find a sensible reasoning for tuning the 
parameters with an accuracy necessary to realize Koide's mass formula. 

We may reverse the argument partially and search for a realistic vacuum within a restricted 
set of configurations. Namely, in view of the high accuracy with which Koide's mass formula is 
realized in Nature, it may make sense to assume that this mass relation is protected by some 
mechanism. (An example of such a mechanism will be given in the next section.) We assume that 
the vacuum configuration satisfies 

($ ) 2 = $ a $ a ; $°GR (34) 

in an appropriate basis allowed by the symmetry. In this case, the relation (28) is satisfied by the 
eigenvalues of <3> [10]; see Fig. 5. Then we minimize the potential of $ within the configurations 
which satisfy this condition. Eq. (34) imposes one condition among the three masses of leptons. 
Apart from the overall normalization of the spectrum, there remains only one free parameter, 
which should be fixed by minimizing the potential. Since the condition (34) or the relation 
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(28) treats the three mass eigenvalues symmetrically, a priori it seems difficult to generate a 
hierarchical spectrum. If we impose the U(3) x £7(3) symmetry to the potential of $, there is no 
vacuum corresponding to a realistic lepton spectrum. We find that in the case of the U(3) x SU(2) 
symmetry, a realistic spectrum follows from the vacuum of a simple potential. 

In the rest of this section we study a classical vacuum of the potential of $ which is invariant 
under the £7 (3) x SU{2) transformation 

$ ^ £7 $ T ; UU ] =0 T = 1. (35) 

Up to dimension 4, there are only 4 independent invariant operators. We parametrize the potential 

as 

V(Q) = Vn($) + V* 2 ($) + V* 3 ($) (36) 

where 

V$i($) = A [tr($ t $) -v 2 ] 2 , (37) 

V* 2 (&) = tr($ f $ $ f $) , (38) 

V$ 3 (§) =e$ 3 tr(<S><S> T <S>*&). (39) 

The 4 independent parameters A, v 2 , e$ 2 and are real. These potentials are classified according 
to the symmetries: Since tr($t $) = | (I)"**!) , V$i is invariant under 50 (18); is invariant under 
SU(3) x S£7(3) x 17(1); V* 3 is invariant under £7(3) x S£7(2) ~ S£7(3) x SU{2) x 17(1). 

We assume the condition (34). When e$ 2 — and A,f 2 ,e$3 > 0, the configuration which 
minimizes U($) under this condition corresponds to a charged lepton spectrum very close to 
the experimentally observed one. Let us describe the details of this configuration. Using the 
transformation (35), any $ can be brought to a form parametrized by 6 real parameters. Without 
loss of generality, we can choose (<E>°, <3> 2 , <3> 3 , $ 4 , $ 6 , $ 8 ) as the real parameters, while $ 1 ,$ 5 ,$ 7 
are set to zero.* Then it is straightforward (but cumbersome) to minimize the potential V$i + V$3 
under the condition (34). One finds the global minimum at the configuration 

«$g,$g)=«b(l, > /T^,xo) 
$ : < (40) 



other $o=0 



where^ 



^ = ( Vi» + 9 )^-(^9-,)V. =ft29gr (41) 



^0 



1 - 3^X0 + 4x 2 g $3 
24 A 



V^l + 0.003656 (e$ 3 /A) ' 



(42) 



There are no other degenerate vacua except those which are connected to $o by the £7(3) x 
SU(2) transformation (35). Note that there is a residual £7(1) T 2 symmetry corresponding to the 
transformation 

$ -> exp(«0T 2 ) $ exp(-i0T 2 ) , (43) 



* Since any $ can be diagonalized by a bi-unitary transformation, $^ = U<&V\ <J> can be brought to an hermite 
matrix by a 17(3) transformation as <&' = V^UQ = V^ d V, i.e. $ Q ' <G R. Noting that ($ 2 ', $ 5 ', $ 7 ') transforms as 
the 3 of the diagonal subgroup SU(2)v C U(3)x 67/(2), we may set <£> 5 ', $ 7 ' = using this transformation. Using 
a residual degree of freedom, which rotates ($ 1 , $ 3 ) as a real doublet of 0(2), we can set = 0. 

^xo is a real solution to the equation 8x 3 + 4x — \J~2 = 0. 
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which keeps the above vacuum invariant. This is a subgroup of U(3) x 577(2). 
The three eigenvalues of $ are given by 

(vi,v 2 , v 3 ) = ^ (VQ + V3x - 3\Jl-xl, V6 - 2\^3x , + V3x + 3^1-x^ 

« v (0.01775, 0.2352, 0.9718). (44) 

The corresponding experimental values read 

(y/m~ e , v^v, s/m~ T ) « (0.01647, 0.2369, 0.9714) , (45) 

where = m e + m /J + m T . We pay particular attention to the value oiv^/vo, which approximates 
the corresponding experimental value with an accuracy of 4x 10~ 4 . Since the constraint (34) treats 
the three eigenvalues symmetrically, some kind of fine tuning should be inherent in this vacuum 
configuration corresponding to a hierarchical spectrum. Indeed, this is reflected to the fact that, 
if the value of v^/vq is varied slightly from the above value under the condition (34), variations 
of vi/vq and V2/V0 are fairly enhanced. (Note that the values of vi/vq and i> 2 /t>o are fixed by 
vz/vq.) As a result, a tiny perturbation to the potential can bring all Vi/vo to be consistent with 
the experimental values. For instance, turning on V$2 with e^/^^ ~ — 6 x 10~ 3 will achieve 
this. This feature is indifferent to details of perturbations: they can be any mixture of V$2, 
higher-dimensional operators, and radiatively induced potentials (log potentials). 

The following comparison may illustrate markedness of the above configuration. When this 
configuration is the zeroth-order vacuum, correct orders of magnitude of m e /m s , m M /m s , m r /m s 
are reproduced if perturbations are sufficiently small. By contrast, when the zeroth-order value 
of V3/V0 is in less accurate agreement with the experimental value, a fine tuning of perturba- 
tive contributions is necessary even to reproduce the mass ratios mj/ms with correct orders of 
magnitude.* 

We find it quite intriguing that the vacuum of such a simple potential, which respects the U(3) 
family symmetry, selects this particular value of v 3 /v very close to the realistic value. Noting the 
relation (27) between the lepton pole masses and the VEV of $ at high energy scales, the above 
feature may suggest that the potential takes a form V(&) V$i($) + V$ 3 ($) at the cut-off scale 
fJL = A. 

At this stage, it is unclear what mechanism protects the condition (34). Furthermore, it is 
unclear why V$ 2 should be so much suppressed compared to V$ 3 , |e$2/^$3| ~ 10~ 2 . Naively, one 
would expect that radiative corrections induce V$2 at least with a similar order of magnitude as 
V$3. In the next section, we will present a possible mechanism or scenario to solve these problems 
(not completely but at least in such a way to circumvent fine tuning of parameters). 

5 A Minimal Potential 

In this section, we present a potential of scalar fields, possibly minimal in its content, which 
realizes (<&) « $ > defined in eq. (40), at its classical vacuum. This is discussed within an effective 
theory which has U(3) x SU(2) family gauge symmetry, valid below the cut-off scale A. The 
assignment to SU(3) x SU{2) x U(l) of the fields, which are already introduced in the previous 
sections, reads 

^l:(3,1,1), e R : (3,1,-1), $ : (3,3,1), y? : (1,1,0), (46) 

tFor instance, if v^/vo ~ 0.9856 (1.5% difference from the experimental value), m e and m M are predicted to be 
the same, m e /ms = m^/ms. 
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with the transformation properties 

if>L^U,p L , e R ^U*e R , $^r/$0 T , p^p, (47) 
U = exp(i6 a T a ) , O = exp(2i6 x T x ) (x = 2,5,7) ; U ' U ] = O T = 1 . (48) 

Furthermore, we assume that above the cut-off scale A there is an SU(9) x U(l) gauge symmetry 
and that this symmetry is spontaneously broken to U (3) x SU (2) below the cut-off scale. 

Let us describe the assignment of the fields to the group SU(9) x U(l). $ is assigned to 
(9,1); its transformation is given by $ a — > £7 Q?/3 $ /3 with a 9-by-9 unitary matrix* U a @ . ipi is 
included in (36, 1) (the 36 is the second-rank antisymmetric representation), which decomposes 
into (6,3,1) © (3,1,1) © (3,5,1) after the symmetry breakdown; similarly e R is included in 
(36, -1). p is a singlet under SU(9) x U(l). 

In order to realize a desirable vacuum configuration, we introduce another field X, which is in 
the representation (45, Qx) (the 45 is the second-rank symmetric representation) and is unitary. 
It can be represented by a 9-by-9 unitary symmetric matrix: 

X a @ — x^ a X ai X^ 1 * = 5 alS ■ X alS U ap X pa U 13 " (49) 

X decomposes into Xg(6, 1, Qx) © X s (6, 5, Qx) © -^^(3, 3, Qx) after the symmetry breakdown. 
See Appendix B for the decomposition of X under U(3) x SU(2). 

We may summarize the essence of how to realize a vacuum, which satisfies eq. (34), as follows. 
If the VEV of X takes a value 

Xf = [diag.(-l, +!,■■-, +1)U = -2 S a0 6f* + 5 a ? , (50) 

an SU(9) x [/(l)-invariant condition 

$ Q X a<i * ¥ = (51) 

reduces to the first condition in eq. (34) at X = Xq. The second condition in eq. (34) can be 
realized by maximizing |$°| upon fixing the value of $ Q *(E> Q and imposing the first condition of 
eq. (34); see Appendix C.l. These conditions can be met at the classical vacuum of the potential 
of the scalar fields under consideration, with an appropriate choice of parameters in the potential. 
We may avoid fine tuning of the parameters, except for the one related to stabilization of the 
electroweak scale. 

In what follows we do not discuss any details of the theory above the scale A. Rather we use 
general properties of SU(9) x U(l) gauge symmetry to infer boundary conditions to be imposed 
at the scale A. We also investigate boundary conditions at this scale required from the low-energy 
side phenomenologically, consistently with symmetry requirements. 

We study the potential and its vacuum of the scalar fields, $, X and tp. First we analyze 
the potential of a specific form (or with a specific choice of parameters of the potential), which 
incorporates an essential part of our model. Later we extend the potential to more general forms. 
The potential we analyze reads 

V($, X) = + V$ 3 + Vxi + V K i + V*xi , (52) 
*For instance, 2tr($t$) = i s invariant under SU(3) x SU(2) x U(l) as well as SU{9) x 17(1). 
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where V$i and V$ 3 are defined in eqs. (37) and (39), respectively, and the other potentials are 
defined by 



V X1 = exi tr(T Q T p T 13 T a ) X afi X pa * , (53) 
V K1 = e Kl \^X^^\\ (54) 
V<s> xl = -ezxi v 2 tr(T a T 13 & T p T a $) X aa * X pp . (55) 

All the parameters of the potential, A, £$3, exi, £k, £<s>xi, v, are taken to be positive. Note that 
since the field X is unitary, it is dimensionless. The physical scale of its VEV is determined by the 
kinetic term of X, which is normalized as f x |(D /t X) Q/3 | 2 . Thus, the physical scale of the VEV of 
X is O(fx)- We choose fx to be much smaller than v (the scale of ($)), such that the spectrum 
of the family gauge bosons is determined predominantly by (<&). (See discussion in Sec. 8.) 
One may verify the following properties of the potential: 

• The global minimum of Vxi is at X = X , defined by eq. (50). Degenerate configurations 
are only those which are connected to X by the 577(3) x SU(3) x U(l) transformation (the 
symmetry transformation of Vxi). 

• Vk\ is minimized if eq. (51) holds. This equation reduces to the condition ($ ) 2 = $ a $ a in 
the case that X = X . 

S0|2 



If X = X , V$xi ~ _ I U P to a term that can be absorbed in V^: 



X=X 



- |<f>°| 2 + — 

8 1 1 18 



(56) 



• If the constraints ($ ) 2 = $ a $ a and $ Q *$ a = 2w 2 (> 0) are imposed, V$xi\x=x is minimized 
when <£> a 's have a common phase 9 £ R, namely e~ t6 <& a £ R for all a; see Appendix C.l. 

• $ = $0, defined by eq. (40), is the classical vacuum of V$3 under the constraints $ a *$ Q = 
2vl ($ ) 2 = $ a $ a and £ R. 

• If the constraints ($ ) 2 = $ a $ a and = 2f 2 are imposed, the first derivative of V$3 
vanishes,''' dV$ 3 /d$ a = <9V$ 3 /<9$ a * = 0, at $ = $ - This is not trivial: In general there 
may be a non-zero derivative in an imaginary direction, since $0 is determined assuming 
$ a £ R. 

• All terms in V($,X) except V* 3 is invariant under S77(3) x 5C/(3) x 17(1) (D C/(3) x 517(2)), 
while the variation of V$ 3 is positive semi-definite at $ = $ . Namely, if $(, = C/i$ t4 
(C/.C/t = U2U \ 2 = 1 ) > v* 3 (^o) > V* 3 ($o); see Appendix C.2. 

Due to these properties, the classical vacuum of V($, X) in the limit £$3, <C £xi is given 
by $ = e^$ an d V = X up to a C/(3) x SU(2) transformation, provided that e^xi/^^3 exceeds 
a critical value to assure the reality condition on $ a (up to a common phase): 

— > 0.02164... . (57) 

tit can be shown, for instance, from the invariance of both V$3 and $0 under the U(1)t^ transformation eq. (43) 
and two Z 2 transformations given by $ Q -> (P. a/3 $^)* with P x = diag.(+l, -1, +1, +1, -1, +1, +1, -1, +1) and 
P 2 = diag.(+l, -1, +1, +1, +1, -1, -1, +1, +1). 
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We note that the definition of v should be modified, including the effect of V$xi, from eq. (42) 
to 



V =V 



53 

144 A 



1 + 



24 



X + 4Xq £$3 



1 + 0.3681(e #xl /A) 
1 + 0.003656 (e* 3 /A) 



(58) 



The operators V$ 3 and V$xi, whose couplings need to be suppressed, are non-invariant under 
SU(9) x U(l). This is a key property of our model which allows us to circumvent fine tuning, as 
we will discuss shortly. As far as the charged lepton masses are concerned, the phase 9 can be 
removed by redefining the phases of tpL and e#. Hence, we set 9 = in the following analysis for 
simplicity.* 

The 0{e^) and 0(e$xi) corrections to the vacuum configuration can be computed. In an 
appropriate basis, these are given by 



5$° = v 



1 3 

+ 



_ 6-5V2x Q 
H =z 



16 e K1 13 sxi J : 78£xi 

5$ 2 = - ^Y^X% 5$°, 5$ 8 = -X 5$°, 



5X 02 = 5X 20 = 2^1 



Xr 



3 5-2^^0 
H 3^ 



13exi 



78exi 



where 



<5X U8 = <5X 8U = 2a;o , 

13exi 

all other <5X Q/3 = 0, 



5 v t v 3 ^ + v 3 ) 
£<s>xi 5 £< J>3 



3 l + 2v^xo-8a;g 

H =7: £<i>xi 



78 5x1 



8 



v, 







8 



-0.06690 e $3 ) 



(59) 
(60) 
(61) 

(62) 
(63) 

(64) 



and ViS are given by eq. (44). Hence, violation of Koide's mass formula is expected to be Oie^jsKi) 
or 0{e<f>/exi)i where Eq> represents e$xi or e$ 3 . The explicit expression of the charged lepton spec- 
trum including the above corrections depends on the precise form of the higher-dimensional opera- 
tor which generates the lepton masses. Naively one expects that e$/£Ki, £&/£xi, £<s>x/£xi ~ 10 -5 
should be satisfied, in order to meet the experimental accuracy of Koide's formula. [Compare with 
the estimates below eq. (94).] 

Next we consider the potential of $ and X in general and examine conditions necessary for 
realizing (<&) $0 and (X) m Xq. Noting that X is unitary, the potential invariant under 
577(9) x U(l) can be written as 



V, 



SU(9)xU(l) 
<5>X 



n,m>0 



C nm ($ a *$ Q ) n \&X^* $' 



I 2m 



(65) 



*The degeneracy of the vacua parametrized by 9 originates from an accidental U(l)<f, global symmetry of the 
potential V($>, X), under which the overall phase of <£> is rotated independently of X . The degeneracy will be lifted 
if we include in the potential operators which break this accidental symmetry. 
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When the coefficients C nm are appropriately chosen (without fine tuning), V^ x can have a 

minimum at > and & X^* $ 7 = 0. These are satisfied by <3> = $o an d X = X . 

All the other operators are non-invariant under SU(9) x U(l). We separate them into three 
categories: those which depend only on X (Vx), those which depend only on $ (V$), and those 
which depend on both $ and X (V$x)- Requirements to each of them are as follows: 

• We can show that the first derivative of Vx vanishes at X = X if CP invariance is preserved; 
see Appendix C.4. This means that, assuming CP invariance, X = X can be the global 
minimum of Vx in a certain domain of the parameter space (spanned by the parameters in 
Vx), without fine tuning of parameters. § 

• Up to dimension 4, V$ consists only of V$ 2 and V$ 3 ; see Sec. 4. Since effects of higher- 
dimensional operators are expected to be suppressed, V$ will be minimized at $ ~ <3>o if 
£$2 <C £$3 <C £k, £x and $ a X^ 13 $^ = 0. Here, Ek/v and Sx/v represent typical magnitudes 
of the second derivatives of V^ x ^ xU ^ and Vx, respectively, at their minima. 

• The contribution of V^ x needs to be suppressed as compared to those of v$ x ^ xU ^ and V X - 
If we can treat V$ x as a perturbation, we may substitute X = X and $ Q *$ Q = Vq in the 
lowest -order approximation. Then V$x becomes dependent only on <3>° and $ x (x = 2, 5, 7). 
The role of the operators dependent only on $° is similar to V$xi, their total contribution 
should not be too small compared to that of V$3 and should enforce the reality condition 
on <3> a ; c.f. eq. (57). The role of the operators dependent on § x is similar to V$ 2 ; in order 
to suppress corrections to the lepton spectrum, contributions of these operators need to be 
suppressed compared to that of V$3. 

Thus, under appropriate conditions, ($) ~ $o can be realized with a more general potential 
than the specific potential eq. (52). Coefficients of certain operators need to be suppressed com- 
pared to the others. One may estimate typical orders of magnitudes of hierarchies required in the 
constraints, from the analysis of the specific potential V(<&,X), which serves as a reference case. 
Moreover, in principle it is straightforward to compute corrections to the vacuum configuration 
similar to eqs. (59)-(63) for a more general potential. 

Let us comment on CP invariance. We may assume that either CP invariance is broken 
explicitly (but weakly) or it is broken spontaneously. In the former case, since there is no observed 
CP asymmetry in the lepton sector, we may assume effects of the explict breaking are very small 
and will not affect our argument given above significantly. In the latter case, since CP asymmetry 
resides only in the Yukawa interaction in the SM, we may attribute the Kobayashi-Maskawa CP 
phase to the VEV of the scalar field, which is presumably existent to give masses to the quarks, 
while keeping all the interactions in the U(3) x SU(2) effective theory CP-invariant. As we do 
not discuss quark sector at all in this paper, this argument is rather ambiguous. In passing, we 
note that all the operators in the potential V(<&, X) [eq. (52)] are CP-invariant; see Appendix C.3 
for the CP transformations. 

Furthermore, the Higgs field tp needs to be incorporated in the potential. Since tp is a singlet 
under U(3) x SU(2), it can be included effectively by replacing the coefficients of the operators 
in the above discussion by functions (polynomials) of tp'tp, e.g. C nm — > C nm (ip<(p). Hence, the 
conditions on the coefficients are the same as above when evaluated at tp'tp = v% w (<^ Vq). On the 

§One example is the case in which Vxi gives a dominant contribution in Vx, although other operators need not 
be suppressed by orders of magnitude. This is because, contributions of other operators cannot create a non-zero 
derivative at X = X , and the position of the global minimum is altered only if their contributions are large enough 
to create a global minimum at another configuration. 
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other hand, the VEV of cp is determined from the same potential after substituting $ ps <£> an d 
X X , whose expansion about the minimum should take a form const. + X v (y?V — ^cw) 2 + ' ' ' ■ 

It is appropriate to regard the conditions discussed above as those to be imposed on the Wilson 
coefficients in the effective potential (in Landau gauge) renormalized at the cut-off scale /i — A. 
Recall that, as we discussed in Sec. 3, we may relate the charged lepton spectrum directly to the 
vacuum configuration of the effective potential at /i — A. The advantage of choosing fi — A is 
that certain fine tuning can be avoided in this way. Let us describe how it works. 

According to the argument above, £$3 is required to be much smaller than Ek or Ex in order 
to suppress corrections to Koide's formula. Furthermore, in Sec. 4 we have seen that e$ 2 should 
be much smaller than £$3 to generate a realistic charged lepton spectrum. Hence, e^/^k and 
£<5>2/£x should be quite small, of order 1CT 5 or less. On the other hand, the 1-loop correction by 
family gauge interaction to the effective potential induces V$2- This indicates that a natural size 
of £$2 is order a\ ~ 10~ 3 or larger within the {7(3) x SU(2) effective theory, assuming the relation 
(26). Thus, in order to realize (<&) ~ $0, a fine tuning of £$2 seems to be requisite (provided 
magnitudes of Ek and Ex are moderate). This argument, however, does not apply at the cut-off 
scale \x = A: Since SU(9) x {7(1) symmetry forbids V$ 2 and V$ 3 in the theory above the scale A, 
both £$2 and e$ 3 are expected to be suppressed at /i — A in the U (3) x SU (2) theory. They are 
determined by the matching conditions at n — A. Radiative corrections within the U (3) x £77(2) 
effective theory essentially do not exist at this scale. 

Another advantage of choosing fx — A in the effective potential is that £77(9) x {7(1) symmetry 
breaking effects on (X) = X are also expected to be suppressed. In other words, the wave 
function renormalizations are common to Xa, X s and X s , to a good approximation. It helps to 
keep the first condition of eq. (34) precise, which follows from eqs. (50) and (51). 

The sizes of Wilson coefficients of operators non-invariant under £{7(9) x {7(1) at n — A depend 
on the dynamics how the breakdown of SU (9) x U (1) gauge symmetry occurs in the theory above 
the scale A. For example, one can imagine cases in which these operators are proportional to 
(powers of) a VEV of some scalar field which breaks £{7(9) x {7(1) symmetry. Then, an operator 
whose dimension is n would have a coefficient of order gA k /M n+k ~ 4 , where g is a combination 
of coupling constants, A is a typical scale of the scalar VEV, k is the power of the VEV, and M 
represents an SU (9) x U (l)-invariant mass scale much larger than A. There are no evident conflicts 
between this naive estimate and the conditions on the Wilson coefficients which we derived above, 
presuming that g can be small but cannot be much larger than unity. For instance, applying the 
estimate to the parameters of V($,X) and V$ 2 , we find 

£77(9) x {7(1) invariant : X,e K1 <0(l), \v 2 < M; (66) 

P A 2 A k f 2 A k 

£77(9) x {7(1) non-invariant : e X \ % £$2, £$3^ j^, £<s>xi ^ j^, (67) 

which are compatible with the desired hierarchy of the parameters e$ 2 <^ £<i>3 ? £<f>xi <^ £ki, £ xi- Of 
course, one should keep in mind that the above estimates are heavily dependent on the dynamics 
above the cut-off scale. 

We may speculate on a possible scenario above the cut-off scale which may lead to (part 
of) the desirable hierarchical relations. Suppose that the symmetry breaking £{7(9) x {7(1) — > 
U (3) x SU (2) is induced by a condensate of a scalar field T?f, which is a 4th-rank tensor under 
£77(9). Indeed if (T$) ~ tr(T a T l3 *T p *T cr ), this symmetry breaking takes place. Through the 
first diagram shown in Fig. 6, the operator e$ 3 tr($ $ T $* $t) may be induced; the double line 
denotes a heavy degree of freedom with an £{7(9) x {7(l)-invariant mass scale M. Since (T?f) ~ 
O(A), the coefficient e$ 3 ~ A/M would be a small parameter provided M ^> A. £$2 is even 
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(T) 




(T a P) 

1-£E± $<*$0$p*$"* ~ £# tr($ $ T $*) 




^ ~ ^- > 0(1) 

X / \ ( T ) 

Figure 6: Speculation on underlying physics that may generate SU(9) x U(l) non-invariant operators. 

more suppressed, since the operator tr(<I>t <3> $t $) cannot be generated by a single insertion 
of (Tpf) at tree level. Either two insertions of (T^f) or a loop correction is necessary, which 
leads to additional suppression factors. The second diagram in Fig. 6 would induce the operator 
sxitr(T a T p T 13 T a ) X al3 X pcr * (together with other operators). Since there is no intermediate 
heavy degree of freedom, the induced coupling Exi, when normalized by A, would be order 1. In 
order to generate e^xi with a desired order of magnitude, we need to suppose a more complicated 
scenario, but we do not pursue this further here, since anyway the argument is quite hand-waving, 
without any explicit model above the cut-off scale. 

To end this section, let us comment on the fine tuning problem in maintaining a large hierarchy 
between the scales, which we mentioned in Sec. 1. In the derivation of the potential of (p, it appears 
unnatural that u ew (<C v) determines the scale, since the natural scales involved in the effective 
potential are A and v before substituting <3> rs $ an d X « X . Currently we do not have any 
reasonable idea on how this hierarchy problem may be resolved. 



6 Inclusion of Another Scalar Field 



Our goal is to generate the spectrum of the charged leptons (m e ,m M ,m r ) such that it satisfies 
Koide's formula with a high accuracy and is proportional to (vf, v\, vf) approximately, where V{'s 
are given in eq. (44). For this purpose, we need to introduce yet another scalar field. This is 
because, if we construct the higher-dimensional operator only from the fields ipL, sr, $, X 
and (p, the corresponding charged lepton mass matrix cannot be brought to a diagonal form with 
any choice of basis allowed by Z7 (3) x 577(2) gauge symmetry. (Note that <E>o is not diagonal.) 
The radiative corrections discussed in Sec. 3 will be altered if the mass matrix cannot be brought 
to a diagonal form, and the QED correction will not be canceled. 

Thus, we introduce a (dimensionless) scalar field £y which is in the (6, 1, Qy) under S77(3) x 
SU{2) x £7(1). It is given as a 3-by-3 symmetric matrix and transforms as Xy — > UT,yU t . 
Consider the potentials 



Vx Y = -£ Y i v a tr (£ y Sy j + e Y2 v 4 tr (£ y £y£ y Sy ) + s Y3 v A 
V*v Y = -e*yitr (£ y <M>t£y$*$ T ) . 



tr ( S y S y 



(68) 
(69) 



We take all the parameters eyi, £y2, £y3, £$yi to be positive. One can show that, for a given $ 
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and in the limit e^yi <C £yi,£y2,£y3, Vt. y + V<s>y. y ^ s minimized at 



• ^ 2 ( j;W 3 r (70) 

Here, U<& is a unitary matrix which diagonalizes $<E>t, i.e., ?7^,$$^?7$ is a diagonal matrix; see 
Appendix D.l. In the case that $ = $ 0) the corresponding unitary matrix is given by 




U = — \ i 1) ; Ul$ Q Uo = ®d 



f V! 

o 2 | . (71) 
\ D 3 ' 



Therefore, we may incorporate £y in the operator to diagonalize the lepton mass matrix. 

As in the previous section, we embed £y in a representation of a larger symmetry group, which 
is valid above the cut-off scale A. We could find a reasonable potential only when we embed Sy to 
a second-rank antisymmetric representation, and this is not possible with SU(9) x U(l). We find 
a way out by enlarging the gauge group. Instead of SU (9) x U(l) we assume that SU (ran) x U(l) 
(n > 4, m > 5) gauge symmetry is exact above the cut-off scale.* Under this symmetry group, $ is 

embedded in the (nm, 1). We denote the field in the latter representation by <E>^ (0 < £ < nm — 1) 
and identify $ ? = for < £ < 8. $ decomposes into a (3, 3, 1) (= $), n - 3 (1, 3, l)'s, m - 3 
(3, 1, l)'s, and (n — 3) (to — 3) singlets after the symmetry is broken down to SU(3) x SU(2) xU(l). 
Similarly X is embedded in the second-rank symmetric representation of SU(nm), denoted by X, 

and X = X^ v for < £, r] < 8. Sy is embedded in the second-rank antisymmetric representation 
of SU(nm), denoted by Y; see Appendix D.2 for the explicit relation between Y and £y. Both 
X and Y are unitary fields. The kinetic terms of X and Y are normalized as f-^ \(D^X)^ V \ 2 and 
fy |(D M F)^| 2 , respectively, where and /y are assumed to be much smaller than v. 

We examine the general potential of $, X and Y which is invariant under SU(nm) x U(l). In 
particular, we would like to see if the potential can be minimized at 

= (o<^<8) (72) 

• \0 (£>8) ' [U) 

X 6n = -2 <5 f0 5 r '° + 6^ , (73) 
Sy = aW T , (74) 

without fine tuning of parameters, where Sy is embedded in Y appropriately. The general poten- 
tial can be written in the following form:''' 

v wgn)*u W= £ c( P „pU^Dn^fe)M^(pD^r- (75) 

Pi,p'i,<ii,q'i>o i=1 

Qtot=0 



*SU(nm) includes SU(n) x SU(m) as a maximal subgroup. Below the cut-off A, the symmetry is broken down 
to SU(3) x SU(2) x U(l), where SJ7(3) and S77(2), respectively, are subgroups of SU(n) and SU(m): SU(S) is 
embedded trivially in SU(n), i.e., the n decomposes into a 3 and n — 3 singlets; SU(2) is a maximal subgroup of 
SU(3)', which is embedded in SU(m) trivially. 

^For instance, the right-hand side of eq. (65), after replacing <!> by $ and X by X, is included in this expression; 
it corresponds to the terms for which q\, <^\,P2iP2iPiiP^i lit = 0- 
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Zi(pi) denote S'f/(nm)-invariant operators* 



zi{pi) = Tr 



t T7\2pi 



pr.F) 



z 2 (p 2 ) = <5>-{X^Y) P2 -<5>*, (76) 

* 3 (p 3 ) = $-(X t -F) 2p3 .X t .$, zM = <b-{Y ] .xf Pi+l -Y ] -<5. (77) 
The summation is constrained to the sector with vanishing U(l) charge by the condition 

Q to t = qi J2 Q^iipi)) -iY. Q( Z M)) = o , (78) 

i i 

where Q(z) represents the U(l) charge of the operator z. Due to complexity of the potential, we 
were unable to clarify if the configuration eqs. (72)-(74) can be a classical vacuum in a sufficiently 
general region of the parameter space spanned by {C(pi,p^, cQ}. We only confirmed this in 
a restricted region of the parameter space: For definiteness, we set (n,m) = (4,5); we consider 
the parameter space spanned by C(pi, p[, for p^p\ < 1 and arbitrary qi,q[, while all other 

C(pi,p'i,qi,q'i) are set equal to zero. In this restricted parameter space, there exists a domain 
with a finite volume (non-zero measure), in which yj^(" m ) xt/ ( 1 ) j s minimized at the configuration 
eqs. (72)-(74) by appropriately choosing Y. Namely, the desired configuration is a vacuum (in 
fact, one of many degenerate vacua) in this domain. See Appendix D.3 for details. This feature 
may indicate that the configuration eqs. (72)-(74) can be a vacuum 

oiv SU_(nm)xU(l) fine 

tuning of the parameters in the potential. 

Operators non-invariant under SU(nm) x U(l) are induced at /i < A. Suppose the following 
operators are induced: 

%,resid = ^ 2 El*T' ( ?9 ) 

V X i, Vm, V^xi as defined in eqs. (53), (39), (55) , (80) 

Vs y , V$s y as defined in eqs. (68), (69) . (81) 

Then, with appropriate hierarchy of the parameters, which we already discussed in this and 
previous sections, we have the configuration eqs. (72)-(74) as a global minimum of the potential. § 
Generally, it depends on the dynamics above the cut-off scale which SU(nm) x [/(l)-breaking 
operators are induced, and a set of operators more general than eqs. (79)-(81) can also lead to 
the same vacuum configuration; see the discussion in the previous section. 
For later convenience, we may take the potential 

V($, X, Y) = V%! + V m + V% ircsid + Vxi + V* 3 + VWi + Vx Y + V^ Y (82) 

with 

^ = \(j^-^j , (83) 

V m = e K1 \&JP' m W\ 2 , (84) 



■tThe dot (•) denotes contraction of SU{nm) indices £, r\, . . . ; X v — X X ■ ■ ■ X , Tr(X) = X^ , etc. 

p 

§There are a number of unwanted massless modes at this minimum. They are included in X, Y and do not 
couple directly to ipL, en, ip, $ and Sy. Although it is straightforward to write down U(3) x S'C/(2)-invariant 
operators which give masses to these massless modes, we do not include those operators, for the sake of simplicity. 
In particular, they do not affect the formulas given in the following discussion. 
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Figure 7: Scalar loop diagram which contributes to charged lepton masses. Dashed line represents all 
the mass eigenstates in $ which couple to the operator of (represented by <8>). 

as a reference potential, instead of V(<&,X) denned in eq. (52). For definiteness, we set (n,m) = 
(4,5). In this case, we obtain the desired vacuum configuration, eqs. (72)-(74), in the limit 
£<s>3,£<s>xi <C £ki,£xi and £<$> Y i <C £yi,£y2,£y3, and with the additional conditions eqs. (57) and 
(144). 

7 Higher— dimensional Operator 

We present a candidate of the higher-dimensional operator which generates the charged 
lepton spectrum. The VEV $o> given by eq. (40), cannot be brought to a diagonal form (in 3-by- 
3 matrix representation) using the U(3) x SU(2) transformation. Hence, the operator such as the 
one in eq. (13) is inappropriate. In terms of the fields which we introduced, a simplest possibility 
may be given by 

Of = ^-^X^X^ Yip e R . (85) 

In this case 2Qx + Qy = is required, such that this operator becomes [/(l)-invariant. Since 
(Xa) ~ — §1 in the basis where (<&) $ [see eq. (105) in Appendix B], the above operator 
can be approximately rendered to the form of eq. (11) by the change of basis, ipi — > UoipL and 
e R — > Uq eR. The corresponding charged lepton mass matrix reads 

25 K^(fj) <JV ew . , 

9V2A 2 ' VM ' 

up to corrections of 0(£^/ek), 0(e$/ex), 0{£^/ £^z) or 0(e$Yi/£Yi), where e$ represents £®xi, 
e$2, £*3) etc. 

We should check whether the radiative correction induced by exchange of $ (Fig. 7) violates 
Koide's mass formula or not. With the reference potential eq. (82), we consider the limit Exi oo 
and £<s>3,£<s>xi,£<i>Yi —* consistently with the assumed hierarchy of the parameters. Physical 
modes of Xa decouple in the former limit. Thus, we consider only and V^ x . One may 
determine the scalar mass eigenstates explicitly and find 




(4tt)< 



l0g <4A^I +1 



-o 



. (87) 
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Since it has a form const, x m i; Koide's formula will not be affected. This result may be more 
non-trivial than one might think at a first glance, since the diagram in Fig. 7 corresponds to 
incorporating the class of (infinite number of) 1-loop diagrams shown in Fig. 2. As a cross check, 
we also computed the coefficient of log/i 2 through renormalization of the operator in the 
symmetric phase (v = 0). 

There are only three physical modes of <3> which gain masses of order v ; these are 5$ a = $ a — <3>q 
which are proportional to $q, $q and iX^ $q. The first mode gives the correction eq. (87), 
while the contributions of the second and third modes cancel. Other modes have masses suppressed 
by e$ 3 , e^xii £$yi, so that their contributions to the loop diagram in Fig. 7 are suppressed. 

In fact the same features apply to the general potential of $ and X, if the parameters of the 
potential satisfies the hierarchical relations required to realize the vacuum configuration $ = $ 
and X = X (as discussed in Sees. 5 and 6). Namely, absence of radiative corrections to Koide's 
formula induced by scalar exchanges can be shown in the limit £$3, e^xi, £$yi, etc. — > (assuming 
that contributions from physical modes of X decouple also in this case). 

As we discussed in Sec. 2 with an example of underlying mechanism, it is assumed that op- 
erators other than Of\ which contribute to the charged lepton masses at higher orders of 1/A, 
are absent (or strongly suppressed) at fi = A. Since these operators are non-invariant under 
SU(nm) x £7(1), sizes of these operators are determined by the physics above the scale A. Within 
the U(3) x SU (2) effective theory starting from this boundary condition, other operators are not 
induced radiatively at lower energy scales and the relation (27) is preserved (see also the discussion 
in Sec. 3).* 

In the limit £$3, <C £ki, £xi and e$yi <C eyi, £Y2, £y3, the root-mass-ratios of the charged 
leptons are given by \fm~ij mo = Vi/vo, where mo = vri\ + m,2 + rriz- They are in reasonable 
agreement with the corresponding experimental values as we have seen in eqs. (44) and (45). 
It would be instructive to see how much corrections are induced to these values by the small 
parameters in the potential. For simplicity, let us compute 0(e$) corrections to the charged 
lepton spectrum (mi,m 2 , m 3 ), corresponding to the potential eq. (82) and the higher-dimensional 
operator eq. (85). The 0(e<s>) corrections read 

0.00209 0.00756\ /0.0312 0.152\ 

e <S>3 + 1 , 




£ki £ xi J \ £ki £xi 

0.00152 0.00406\ /0.0227 0.0833\ 

)s m + + e^xi , (89) 

£ki e xi J \ £ki £xi / 

0.000406 0.00112\ / 0.00607 0.0229\ . . 

+ U*3 + • (90) 

£ki £xi J \ £ki e xi J 

As can be seen, the magnitude of the correction is larger for smaller mass eigenvalues, reflecting 
the nature of a hierarchical spectrum, as we discussed in Sec. 4. Comparing to eqs. (44) and 
(45), one finds constraints on typical orders of magnitude of the parameters as e^/eki P$ 10°, 
^$3/^x1 ~ 10~\ £$xi/£ki ~ 10 _1 , e<s>xi/£xi ~ 10 -2 , provided there is no correlation or fine tun- 
ing among these parameters, or with 0(e$Yi/ 'sn) corrections. On the other hand, the overall 
normalization, 

m 1 + m 2 + m 3 = ^0 i 1 + corr -J > ( 91 ) 

*A simpler operator such as ipL $ ^^YfeR would be inappropriate for a candidate of even though it gives 
the desired spectrum at tree level: This operator induces a mass matrix 5Mi oc 1 radiatively, upon contraction of 
$ and & . 
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is subject to radiative corrections induced by electroweak gauge interaction (including QED), fam- 
ily gauge interaction and scalar exchanges, in addition to the 0{e§) and 0(e^Yi/^Yi) corrections. 

We define the following quantity as a measure of the degree of violation of Koide's mass 
relation: 

A = -^-^ — _ i . 92 

3 (mi + m 2 + m 3 ) 

This quantity vanishes if Koide's relation is satisfied. With the reference potential and Of \ the 
0{e<z>) correction reads 

ZZ-lby/2x \_ 67- 75V2X0 
H tz H — £<s>xi (93) 



%Eki 65exi / 390 exi 

/ 0.00523 0.0171 \ /0.0781 0.346 \ 

~ U*3 + + £$xi ■ (94) 

V £ ki £xi / \ £ki £xi / 

Comparing to the present experimental value A exp = (1.1 ± 1.4) x 10~ 5 , we obtain constraints on 
typical sizes of the parameters more stringent than the previous ones: e^z/ e Ki ~ 10~ 3 , e^/exi ~ 10~ 3 , 
£<s>xi/£ki ~ 10~ 4 , e<pxi/^xi ~ 10~ 4 . 

It is easy to adjust the root-mass-ratios a/ mj/m to be consistent with the current exper- 
imental values without violating Koide's relation, as we discussed in Sec. 4. For instance, it is 
achieved by incorporating V$ 2 with £$2/s$3 ~ — 6 x 10~ 3 into the potential. 



8 Relevant Scales and Further Assumptions 

Let us discuss the energy scales, A, i>3(~ v ), fx, involved in the present model. It would be 
unnatural if there is a large hierarchy between v 3 and A, or between f x and v 3 . As we speculated 
in Sec. 3, the scale of U(3) symmetry breaking, typically given by t> 3 , may be at 10 2 -10 3 TeV, 
such that the QED correction is cancelled within a scenario of unification of the electroweak 
SU(2)l and family U(3) symmetries. There are two indications that the cut-off scale A and the 
U(3) symmetry breaking scale v% are not too far apart. One indication is the importance of the 
universality of the U(l) and SU(3) gauge coupling constants eq. (14). This universality may be 
protected above the cut-off scale by embedding 577(3) and £7(1) in a simple group. If v 3 is very 
different from A, however, the values of the two coupling constants at // ~ 1*3 would become too 
different. Another indication consists in the relation (86), from which one derives 



v 1= ( 9V2m T \ V2 _ 1 
A [25K^av ew ) ~ i7v^W^ 



(95) 



up to electroweak corrections, etc. a is smaller than 1/2 and is expected to be not very much 
smaller. Although it depends on the mechanism how the higher-dimensional operator is 
generated, if is not large (as one naively expects), hierarchy between v$ and A is mild. As 
for the scale fx, it is required to be smaller than ($) in order not to alter the spectrum of the 
family gauge bosons. Numerically f x ^ 0.3 V\ ~ 0.005 v 3 is required from the present constraint 
on Koide's formula. 

There are a few more assumptions implicit in the present model, which we have not discussed 
so far. We assume that the SU (2) family gauge symmetry is broken spontaneously at a scale higher 
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than (<&). This is required to protect the symmetry breaking pattern eq. (22), which constrains 
the form of the radiative correction by the Z7" (3) gauge bosons. To achieve this, we need additional 
fields or dynamics, such as an SU(2) doublet scalar field whose VEV breaks 577(2). As yet, we 
have not succeeded to incorporate such a mechanism consistently into our model. Here, we simply 
assume that the breakdown of SU (2) has occurred without affecting the properties of our model 
described above. 

We also assume cancellation of gauge anomalies and decoupling of unwanted fermions. Namely, 
we assume cancellation of anomalies introduced by the couplings of fermions to family gauge 
bosons, at the scale where U (3) x 577(2) symmetry is unbroken. This means that we need fermions 
other than the SM fermions. Fermions other than the SM fermions are requisite in our model 
also because tpL and e# are embedded into larger multiplets of SU{nm). At lower energy scales, 
all the additional fermions are assumed to aquire masses of order v or larger, so that they 
decouple from the SM sector at and below the electroweak scale. Only the SM fermions remain 
at these scales. Presently we do not have a model which fully explains these features. 



9 Lepton Flavor Violating Processes and Other Predic- 
tions 



A most characteristic prediction of the present model is the existence of lepton-flavor violating 
processes induced by the family gauge interaction. In the scenario, in which the U (3) family gauge 
symmetry and SU(2) L weak gauge symmetry are unified at 10 2 -10 3 TeV scale, the family gauge 
bosons have masses of the order of the unification scale. 

As it is clear from eq. (21), flavor violating decays of a charged lepton with only charged leptons 
and/or photons in the final state, such as /i — > 3 e or /i — > ej, are forbidden. Flavor violating 
leptonic decays which involve neutrinos, such as fx~ — > e"v e v il , are allowed, but the present 
experimental sensitivities are very low. Presumably, the most sensitive process is K L — > e /i, 
although we need to make assumptions on the quark sector. For instance, assuming that the 
down-type quarks are in the same representation of U(3) as the charged leptons, and that the 
mass matrices of the charged leptons and down-type quarks are simultaneously diagonalized in 
an appropriate basis, this process is induced by an effective 4-Fermi interaction connecting the 
first and second generations: 



-(1,2) 1 

4/ 2 K + vl) 



(d Yl5 s + s 7^75 d) (e 7^75 fx + fx 7^75 e) 



(d Y s - s Y d) (e 7„ fj, - \i 7„ e) 



+ 



(96) 



We find 



T{K L -> en) 



™>l™>K L fK 
167Tf| 



(97) 



Comparing to the present experimental bound Br(K L — > e/i) < 4.7 x 10 [2], we obtain a limit 
V2 ~ 5 x 10 2 TeV. Naively this limit may already be marginally in conflict with the estimated 
unification scale in the above scenario. We should note, however, that this depends rather heavily 
on our assumptions on the quark sector. In the case that there exist additional factors in the quark 
sector which suppress the decay width by a few orders of magnitude, we may expect a signal for 
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K L — > e/i not far beyond the present experimental reach. Similarly the process K + — > 7r + e~/i + 
may also be observable in the future. 

Another interesting observation, although it is much more model dependent, is the following. 
In order to stabilize Koide's formula, in our model, it is necessary to suppress SU(nm) x U(l) 
non-invariant operators in the potential of $. This indicates that $ includes physical modes which 
are much lighter than v ~ 10 2 -10 3 TeV. In particular, the lightest one, being singlet under the SM 
gauge group, may decay into leptons through the family gauge interaction or the operator 
with a significant branching ratio. Hence, if this lightest scalar boson happens to be produced at 
the LHC, an excess in multi-lepton final states may be observed. 

10 Summary and Discussion 

In this paper, we propose a model of charged lepton sector, in the context of an EFT valid below 
the cut-off scale A, which predicts a charged lepton spectrum consistently with the experimental 
values. In particular, we implement specific mechanisms into the model, such that the spectrum 
satisfies Koide's mass formula within the present experimental accuracy. In this model radia- 
tive corrections as well as other corrections to Koide's formula are kept under control, and this 
feature primarily differentiates the present model from the other models in the literature which 
predict Koide's formula. By studying within EFT, we circumvent many problems, at the price 
of introducing the cut-off scale at 10 2 -10 3 TeV scale, while non-trivial relations between family 
symmetries and observed charged lepton spectrum can still be investigated. 

In our model, we adopt a mechanism, through which the charged lepton mass matrix becomes 
proportional to the square of the VEV of a scalar field <3> [10]. On the basis of this mechanism, 
we incorporate two new mechanisms in the model which are worth emphasizing: 

(i) The radiative correction to Koide's formula induced by family gauge interaction has the same 
form as the QED correction with opposite sign. This form is determined by the symmetry 
breaking pattern eq. (22) and the representations of ipL and cr. Within a unification scenario, 
cancellation of the QED correction can take place. 

(ii) A charged lepton spectrum, which has a hierarchical structure and approximates the experi- 
mental values, follows from a simple potential V$3, under the condition that Koide's formula 
is protected. 

Existence of such simple mechanisms may indicate relevance of U(3) x SU(2) family gauge sym- 
metry in relation to the charged lepton spectrum. 

Our model is constructed as an effective theory valid below the cut-off scale A respecting 
this symmetry. We introduce scalar fields $, X and Y as multiplets of SU(nm) x U(l), in 
which U(3) x SU(2) is embedded. It is assumed that SU(nm) x U(l) is spontaneously broken to 
U (3) x SU (2) below the scale A. We minimize the potential of the scalar fields and determine its 
classical vacuum. The charged lepton masses are related to the VEVs of the scalar fields at scale 
jj, — A, mf olc oc fj(A) 2 ; at this scale radiative corrections to the VEVs essentially vanish within 
the effective theory. Then, the mass matrix of the charged leptons are given in terms of the VEVs, 
such that Koide's mass formula is stabilized, and that the spectrum agrees with the experimental 
values. This is achieved formally without fine tuning of parameters in the model, except for (a) 
the tuning required for stabilization of the electroweak scale v ew , and (b) the tuning required for 
the cancellation of the QED correction, that is, realizing o:f = \a at relevant scales. We argue 
that the latter tuning can be replaced by a tuning of the unification scale, within a scenario in 
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which U(3) family gauge symmetry and SU(2) L weak gauge symmetry are unified at 10 2 -10 3 TeV 
scale. 

In addition our model may contain following fine tuning. We were unable to explore the pa- 
rameter space of the SU(nm) x [/(l)-invariant potential y^^^ rn ) xC/ ( 1 ) sufficiently, due to technical 
complexity. It may be the case that certain fine tuning is necessary to realize the configuration 
eqs. (72)-(74) as a classical vacuum. 

Evidently the present model is incomplete, since it is restricted to the charged lepton sector. 
The model should be implemented in a larger framework which incorporates at least the following 
aspects missing in the present model: (i) Including the quarks and explaining the masses and 
mixings of the quarks and neutrinos; (ii) Cancellation of anomalies introduced by the couplings 
of fermions to family gauge bosons; (iii) Unification of U(3) and SU(2) L gauge symmetries at 
10 2 10 3 TeV scale. Possibly these problems are solved simultaneously in some model, and one 
anticipates that such a model would necessarily contain a large number of new particles, for the 
following reasons: (a) all the particles are embedded into multiplets of large groups, especially 
if one also requires to unify hypercharge U(l) and color SU(3) gauge groups together with U(3) 
and SU(2) L ; (b) additional fermions are necessary to cancel anomalies; and (c) scalar fields would 
be necessary to give masses of order ($) to fermions (apart from the SM fermions) through their 
VEVs [19]. 

Although our model predicts a realistic lepton spectrum, in fact many of the questions are 
simply reassigned to physics above the cut-off scale and remain unanswered: While we replaced 
the conditions on the lepton spectrum by the boundary conditions of the effective potential, we 
do not address which dynamics leads to these boundary conditions. (Only a speculation is given.) 
We may nevertheless state that not only did we circumvent fine tuning but also the problems 
actually simplified. The required boundary conditions are certain hierarchical structure among 
the couplings of the effective potential. These conditions would be simpler to realize than, for 
instance, to realize Koide's relation among the lepton Yukawa couplings with 1CT 5 accuracy a 
priori 

Phenomenologically our model predicts existence of lepton violating processes at 10 2 -10 3 TeV 
scale, assuming the unification scenario at this scale. The processes K L — * fie and K + — > n + e~ fi + 
are expected to be sensitive to the predictions of our model, although we need additional assump- 
tions on the quark sector. Stability of Koide's formula indicates existence of light modes in $, and 
the lightest mode may decay into leptons with a significant branching ratio; they may generate 
an interesting signal at the LHC. 

It is unlikely that the present model describes Nature correctly to the details, since we can 
easily construct variants of the present model with similar complexity. Overall, the present model 
is rather complicated, and the source of complexity is conspiracy to realize Koide's formula with 
a high accuracy. Hence, we place more emphasis on the major mechanisms incorporated in the 
model, which look appealing and may reflect physics that governs the spectrum of the charged 
leptons. 
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Appendices 



A Generators of U(3) 

The generators for the representation (3, 1) of U(3) ~ SU(3) x U(l) are given by 



T u = -= I 1 I , T 1 = - I 1 I , T 2 



T 3 = - -1 , T 4 = - , T 5 = - , (98) 





f = - 1, T = - —% , T 8 = 




which satisfy eq. (9). 



B Decomposition of X under 17(3) X SU(2) 

X, which is in the (45, Q x ) of 5C/(9) x U(l), decomposes into X%(6,1,Q X ) © Xf (6, 5, Q x ) © 
X^(3,3,<5x) under SU(3) x SU(2) x C/(l). Explicitly they can be constructed as follows: 

X^^X^T^T^ °& X°* = 4X ikijl T? i T l l, (99) 
X s \ = 1 

1 



. , \X ik .ji ± X ki .ji ± X ik; ij + X ki; ij) (ioo) 



(ioi) 



(Xa) Cnjl(XA)ik;jl , (102) 

(Xg)ik = (Xs)ik;mm , (103) 

(Xg)ik;jl — (Xs)ik;jl — g (^slikimm • (104) 

When (X) = X , the corresponding VEVs of Xa, Xg and Xg are given, respectively, by 

5 

(X A ) mn — — - 5 mn , (105) 

(^) ifc = ^(J ifc , (106) 

= Y^Wkj + SijSki) - j^Sikdji ■ (107) 



31 




Figure 8: Eq. (110) corresponds to a nonagon with a fixed length of circumference in the complex plane. 
If we maximize \z°\, the nonagon collapses to a line; after overall phase rotation, all z a, s can be made 
real, where z° = v q > and z a < 0. 

C Properties of X) 

C.l Maximizing |<3>°| 2 

We impose the conditions 

= 2vl >0, ($ ) 2 = (108) 

If we maximize |$°| 2 under these conditions, all $ Q 's can be made real simultaneously by a common 
phase rotation. Namely, there exists a phase 9 such that e~* 9 $ Q £ R for all a. 
Proof. Let 

2 ° = ($ ) 2 , ^ = -($ a ) 2 . (109) 

Then z a, s satisfy 

8 8 

These equations represent a nonagon with a fixed length of circumference in the complex plane. 
\z°\ = |$°| 2 is maximized when the nonagon collapses to a line, where all z a, s are parallel to one 
another and antiparallel to z° in the complex plane with 



\z°\ 



a=l 



See Fig. 8. 



C.2 Variation of V^ 3 at <3> = <3> 

The variation of V$3, defined by eq. (39), is positive semi-definite at $ = $o under SU(3)xSU(3) x 
17(1) transformation. Namely, if $' = C/i$ f4 with u i u i = W| = !> ^W^o) > V* 3 ($ )- 
Proof: Let 

$1 = Wj^o = Mot4 ■ (H2) 
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Then, noting 



($°) 2 = = Vo 2 , (113) 

$1 satisfies the same relation: 

($0) 2 = = v 2 . (114) 

Since $q 's are real, $"'s are also real. According to Sec. 4, $ is a configuration which minimizes 
V$ 3 (<3>) under the condition ($ ) 2 = $ a <l> a = v$ and $ a G R. Therefore, V$ 3 ($i) > V$ 3 (<E>o). 
Since <3>i and $q are connected by a £7(3) x SU(2) transformation, V$ 3 ($i) = V$ 3 ($' ). It follows 

Due to this property, Vxi + Vk + V$ 3 is minimized at X = X and $ = <£> under the constraint 
fGR and = 2w 2 . 



C.3 CP transformations of <J> and X 

CP transformation of $ is defined by 

(CP) $(x) (CP) 1 " = $(Vx)* , (115) 
where Pa; = (a; , — x). Equivalently, 

(CP) <5> a (x) (CP) f = [C af3 ^(Vx)}* (116) 

with 

C a(B = [diag.(+l, +1, -1, +1, +1, -1, +1, -1, +l)] Q/3 . (117) 
Similarly CP transformation of X is defined by 

(CP) X a ?(x) (CP) 1 " = [C aa ' X a 'P'(Vx) C w ]* , (118) 

or 

(CP)(X A ) l] (CPy = (X A )* J , (119) 
(CP) (Xg)ij (CP) 1 " = (Xl)*j , (120) 
(CP) (Xg) ik .ji (CP) 1 " = (Xg)* k;jl . (121) 

CP transformations of other fields are the same as those of the SM. 

For example, Vxi defined by eq. (53) is CP-invariant. An example of CP non-invariant 
operator is 

v =^(f-n ( 122 ) 

with 

/ = [tijkUmn (X A )u (X A )jm (X A )kn\ [^i'j'k'^l'm'n' (^s)i'J' (^s)j'm' (Xg)k'n']* ■ (123) 
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C.4 Stability of V x at X = X 

When Vx(X) is invariant under U(3) x SU(2) and CP, its first derivative vanishes dVx/dX af3 

OV x /dX a/3 * = at X = X . 

Proof: 



SV; 



X 



dV x /cv , , SV^ rAy dV x 



+ 



d{X A ) i3 
dV x 



(SX. 



A)ij 



(sx, 



1 \* 
SJij 



+ 



dV x 



(SXg) ik .ji + 



(6Xl)ij 
dV x 



5n* 

S7ifc;j7 • 



(124) 



Due to the residual (2)y symmetry of the VEV (X) = X , the differential coefficients evaluated 
at X = X Q take following forms: 



x 



x 



x 



dV; 



X 



W ' d{X A % ' d{Xl) i3 ' 9(Xi 



dV x 



OC (L- , 



d(X b s ) ik . 3 i ' <9(Xf )* fc;j7 



(125) 
(126) 



X=X 



where Cj's are constants. Substituting to eq. (124), we have 



SV- 



x 



C[ (SX A ) U + C 2 {5X A )l + C 3 (SXDh + C' A {5Xl% 



x=x 



+ C 5 {SX%) ik , ik + C 6 (5Xlr ik]tk + C 7 (5Xl) tk , kt + C 8 (5Xl)* k;kl , (127) 
where we used (SXg)ajj = 0. 



An arbitrary infinitesimal variation of X, which is symmetric and unitary, can be parametrized 



by 



with 



x a(S + sx a(S = w aa ' x a 'P' ww 



w 



( l + ie§ i(e 



lie 



08 



ie 



08/ 



«(4 + ie is) 



(128) 



(129) 



\ ^( e 08 ^ e 08) "i(. e 18 ^ e 18) '■' 1 + ^ e 88 / 

neglecting 0(e 2 ) terms. An explicit calculation shows that, for a variation X = X + 5X, (5X A ) ii: 
(SXg)a, (SXg)ik-ik and (SXg) ik - ki depend only on e^ a for < a < 8. (In this proof, no sum is 
taken over a in e^ a without explicit summation symbol ^2 a -) Hence, 



SV; 



x 



M 



d 



X=X a=0 Ut aa 



V X (X + 5X) 



(130) 



e « , e < „=0 

a/3' ap 



On the other hand, applying CP transformation eq. (118) to X = Xq + SX, one finds that Xo is 
CP-even, whereas all the coefficients of ef a in SX are CP-odd. This means, if Vx is CP-invariant, 



d 



■V X (X + SX) 



= for < a < 8 , 



(131) 



so that the first derivative vanishes, SV x \x=x = 0. 



34 



D Y, E y and Their Potential 
D.l Minimum of Vz Y + V$s v 

We show that Vs Y + V$s y , given by eqs. (68) and (69), is minimized at the configuration eq. (70) 
in the limit e^yi <C £yi,£y2,£y3- 

It is known [18] that, for £yi,£y2,£y3 > 0, V^ Y is minimized at 



*r = °VV- : ^,/ 2(£y J' 3£r3) . (132) 



where C/ is an arbitrary 3-by-3 unitary matrix. Let 

/ u\ \ 

C/J$$ t C/$=[ M2 \=Ml ; Mi >0. (133) 
V u§ / 

We assume that all u^s are different. Substituting eqs. (132), (133) to V$£ y , it is expressed as 

V*z Y = -e* Y i<J 2 tr(W* M 2 WM 2 ) ; W = u\, U U T U% . (134) 

W is unitary. Define 

M 2 d =A a T a , W^M 2 d W = B a T a ; A a , B a G R. (135) 

Then = £ a £ a , since tr \(W* M 2 d wf] = tr(Afj). Hence, V* Ey = -± e*yi aM Q i3 a is 

minimized when „4. a = This means W = Ud and L^ 7 " = U$UdU$, where Ud is an arbitrary 
diagonal unitary matrix defined in eq. (24); it can be absorbed into a redefinition of U& as U'^ = 

D.2 Relation between Y and Ey 

Y is in the (nmCiiQY) under SU{nm) x £/(l), where nrn Ci stands for the second-rank antisym- 
metric representation of SU{nm). Y is defined to be unitary. Thus, 

y*» = -y* ; Y (v Y (ri * =5^. (136) 

The indices take values < £, 77, £, • • • < ?im — 1. 

An orthonormal basis of n-by-m matrices is denoted by {T^} with the normalization condition 

tr(T e V) = tr^T^) = 1^. (137) 
In particular, the first 9 bases are taken as 



4 \ \f l,J - 3 (o<e<8). ass) 

otherwise 
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We may identify Ey(6, 1, Q Y ) embedded in Y as follows. 



Y SA \ I 
Y AS 



= | (Yik-ji ± Yia-ji T Y ik;lj - Y kt ,^ (140) 

ik;jl 

= \(Y tk . s ±Y kv ^ , (141) 

(E Y )ik = {YsA)ik- A5 for 1 < i, k < 3 . (142) 



D.3 A vacuum of the SU(nm) X U(l)— invariant potential 



We analyze a vacuum configuration of the SU(nm) x [/(l)-invariant potential y^^™ m ) xC/ ( 1 ) gj ven 



by eq. (75). We restrict our analysis to the case (n,m) = (4,5) and consider only C(pi,p' i ,q i ,q' i ) 
for Pi, < 1 and arbitrary qi,q[, while all other C(pi,p' i ,qi,q' i ) are set equal to zero. In this 
restricted parameter space spanned by {C(pj,p^, <^)}, we examine if the configuration given by 
eqs. (72)-(74) can minimize y^(^_ m ) xC/ ( 1 )_ We assume that the £7(1) charge vanishes, 

Q tot = Qi J2 Q(zi(Pi)) -q' t J2 QWJ) = > ( 143 ) 

only in the sector for which X]j(<?« + ?D > 1- This is not a strong condition: Except when and 
satisfy specific relations, this condition is met. 

We have checked the following two properties. (I) At each point of the parameter space, the 
first derivative of y^(^ m ) xC7 ( 1 ) vanishes at the configuration eqs. (72)-(74), if 

a<\. (144) 

(II) The configuration eqs. (72)-(74) minimizes V^ l ^y n ^ xU ^ 1 \ if the condition (144) is satisfied 
and at each point in a hypersurface S in the parameter space; the hypersurface S is defined 
by the condition C(pi,p' i ,qi,q' i ) > if pi = p\ and qi = q[ for all i, while C(pi,p' i , qi, q[) = if 
Pi Pi or Qi f° r an y i- These two properties (I) (II) ensure that, of each point in S, there 
exists a neighborhood, which has a non-zero volume, and in which y/£^(™' Tt ) x ' 7 ( 1 ) j g minimized 
by the configuration in question. Namely, there exists a finite volume (non-zero measure) in the 
parameter space (at least) in a neighborhood of S, in which the desired configuration becomes a 
vacuum. 

The above properties (I) (II) are verified in the following manner. It suffices to show that all 
Zi{pi) for pi < 1 can be brought to zero simultaneously at the configuration eqs. (72)-(74) by 
appropriately adjusting components of Y except for Ey. In fact, in this case, y/£^_ m ) xC/ ( 1 ) as 
well as its first derivative vanish at any point of the parameter space. Thus, property (I) follows. 
Since y^^^ m ) xC/ ( 1 ) > o in S*, the property (II) follows as well. We have checked numerically that 
all Ziipi) can be brought to zero at the configuration eqs. (72)-(74) by explicitly constructing the 
corresponding Y for a given value of a. This turned out to be possible (at least) if the condition 
(144) is met, since there are quite large degrees of freedom in the choice of Y . (If a is too large, 
it conflicts the unitarity condition of Y .) 
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